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Abstract 

We discuss the AdS / GET duality from the perspective of integrable systems and estabhsh a direct relationship 
between the dimension of single trace local operators composed of two types of scalar fields in A/" = 4 super Yang- 
Mills and the energy of their dual semiclassical string states in AdS^ x S^. The anomalous dimensions can be 
computed using a set of Bethe equations, which for "long" operators reduces to a Riemann-Hilbert problem. We 
develop a unified approach to the long wavelength Bethe equations, the classical ferromagnet and the classical 
string solutions in the SU (2) sector and present a general solution, governed by complex curves endowed with 
meromorphic differentials with integer periods. Using this solution we compute the anomalous dimensions of 
these long operators up to two loops and demonstrate that they agree with string-theory predictions. 
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1 Introduction 



Gauge/string duality is an old and fascinating subject that appears in a variety of situations 1 . This subject 
encompasses a large circle of ideas, the most promising of which, perhaps, is the relationship between strings 
and planar diagrams in the large- limit 7. A well known example of such a duality is the matrix model 
description of two-dimensional quantum gravity and noncritical strings [SJEl- Matrix models 5,, together with 
two-dimensional QCD [S], are two instructive examples in which the large- iV limit is solvable. 

A much more complicated case of this large- duality is the AdS / CFT correspondence, an asserted equiva- 
lence of a four-dimensional gauge theory with J\f = A supersymmetry and type IIB string theory on the product 
space AdS^ x 5*5. Among other predictions, the AdS/CFT conjecture relates the dimensions of gauge- invariant 
operators with the energies of particular closed string states propagating in the AdS^ x background [7| El EI ■ 
At present, it has only been possible to test this conjecture for a limited class of operators. The reason, of 
course, is that the string calculations are trustworthy for large 't Hooft coupling A — g'^N, while the gauge 
theory calculations are reliable if A is small. 

There are two basic sets of operators which get around this mismatch of the coupling strength. The first 
set contains the chiral primary operators and their descendants. Here one may rely on nonrenormalization 
theorems to test directly the correspondence |7| |H1 E] ■ The second set consists of operators with large global 
charges, which are dual to semiclassical string states. An example of this operator type are the so called BMN 
operators 1101. Here, one starts with a chiral primary made up of a large number of complex scalar fields 
(Tr#^ for i ^ 1), adds to the operator a handful of "impurities", which can either be other scalar fields or 
covariant derivatives acting on the $ fields, and then sums over the positions of the impurities weighted by 
some phases. On a superficial level, the dual closed string can be visualized as a chain of scalar fields inside a 
trace, where the impurities on the chain correspond to excitations of the string. The precise AdS dual of the 
chiral primary is a point-like string orbiting a geodesic of with angular momentum L, large enough for the 
classical approximation to be accurate The nearby spacetime of the string trajectory has a plane wave 

geometry in which the world-sheet string action in the light-cone gauge can be formulated as an infinite set of 
massive oscillators '121 113| . Therefore a direct comparison can be made between the anomalous dimensions of 
the perturbative gauge theory |1U[I14|[T^ and the energies of the string oscillations in the plane wave. 

Since the BMN operators are assumed to have a small number of impurities, the strings dual to these 
operators are almost pointlike. Once the string has significant stretching, the exact string spectrum is unknown. 
Nevertheless, following UJ, a large number of classical string solutions in AdS 5 x have been constructed 
[IBlinilIlliniliniEIll22E3ElE3- These string configurations are described by solitons of the world-sheet 
sigma-model, all of which have large actions and therefore correspond to semiclassical states that are high in the 
energy spectrum of the string'^ . A potential difficulty in comparing string solitons to the supersymmetric Yang- 
Mills (SYM) operators is that the operators are necessarily large, i. e. they must contain many constituent fields 
in order to carry large quantum numbers. The computation of the anomalous dimension for large operators, 
even at one loop, is hindered by operator mixing which lifts the degeneracy of the classical scaling dimensions. 
The degeneracy grows exponentially with the size of the operators and the mixing becomes more and more 
complicated. Fortunately, the mixing matrix can be identified with the Hamiltonian of an integrable spin 
chain 23 |2H1 HHj- One can then use powerful techniques of the algebraic Bethe ansatz |SI3] (see the up to 
date formulation of this method in |Hlj ) to diagonalize the mixing matrix and to compute the anomalous 
dimensions. Explicit calculations have led to a remarkable agreement with string theory predictions in several 
cases ISSESIEIEI]. 

comprehensive review of classical string solutions in AdSs X can be found in I2til . 
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Much of the recent progress focused on single trace operators composed of two types of complex scalar fields, 
$1 and $2, of the form 

Tr (*l$l$l$2*2*l4'2*2*l*l*l*2 ■ ■ ) ■ (11) 

It was shown in 27 that the one loop dilatation generator when acting on operators of this type is equivalent 
to the Hamiltonian of the XXX Heisenberg spin chain. Hence, we will call operators with the form in 
XXX operators. These operators are naturally identified with the states of the spin chain by associating, say, 
$1 with up spins and #2 with down spins. For instance, the operator is mapped to the state 

iTTTiiTiiTTTi-..) e (c^)^"^. 

The one loop dilatation operator of A/" = 4 planar SYM theory was found in , and its two loop correction 
was proposed in ,38j. It has the following form in this basis 

A ^ 

(1.2) 

((1 - ■ ^1+2) - 4 (1 - • CTl + l)) + 0(A=^), 

1=1 

where ct; are the Pauli matrices which act on the spin at the l^^ site of the chain. The dilatation operator is 
not diagonal in the basis of monomials and mixes simple operators of type Hl.l|l . It can be diagonalized, at one 
loop (first line in (|1.2|l ) by using the Bethe ansatz technique 30 . 

On the string side, restriction to XXX operators corresponds to a string moving on the subspace C 
and localized to the center of AdS^. In \6b\ it was shown how to directly compare the equations of motion for 
a classical spin chain with the string sigma model equations of motion by taking a "nonrelativistic" limit of the 
effective coupling. The classical chain naturally corresponds to the long wave-length limit of the quantum chain 
|36l IS7] , which is precisely what was considered in |32[ I33L . 

An obvious question is how the integrability of the spin chain relates to integrability of the sigma model. 
There is a body of evidence that the integrability of the mixing matrix for the XXX operators is not an 
accidental symmetry of the one loop approximation. The XXX operators mix only with themselves to all orders 
in pertubation theory and in this subsector the integrability definitely extends to two loops f3S| , to three loops 
[SHI 1^ ilO., M4, and possibly to higher loops 021 ■ On top of this, it has been verified that the superstring 
sigma-model in AdS^ x is classically integrable |43l I44j . Furthermore, it was shown how solutions on the 
sigma model reduce to a classically integrable mechanics, allowing one to explicitly solve for the string motion 
The integrable structures on the two sides of the AdS/CFT correspondence should be related in some 
fashion. On an algebraic level, the relation among the integrals of motion was discussed in 45 . Subsequently, 
Arutyunov and Staudacher were able to make a direct comparison between certain string solitons and large 
XXX operators in SYM, showing that the entire integrable hierarchy matches in the small coupling limit |46) . 
This matching of hierachies has also been extended to certain classical solutions of the 5[/(3) chain j47| . 

The main goal of this paper is to develop a unified approach to the Bethe ansatz, string solitons and the weak- 
coupling "non-relativistic" reduction of the sigma-model. We will discuss in detail mainly the XXX subsector, 
though we will also mention a set of string solutions that are not dual to this class of operators. 

In comparing perturbative gauge theory results to string theory predictions, it is important to note that gauge 
theory and string theory compute different expansions for the anomalous dimension. The string calculations are 
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accurate when both the bare dimension of the operator, L, and the 't Hooft couphng A are large, but the ratio 
A/L^ is finite. Since the relevant parameter is A/L^, it is possible to consider strong coupling on the string side 
and still have A/L^ small. One might then expect to make a direct comparison to perturbative gauge theory. 

However, it was argued by Serban and Staudacher 42 that nonpcrturbative terms could complicate the 
situation. These authors showed that a particular long range spin chain would reproduce the two-loop and 
three-loop dilatation operators of j38i .SQt .40) if one identifies X/L^ with a parameter of the spin chain. They 
then went on to show that the operators dual to a folded and a circular string match at two- loops, but fail to 
do so at three-loops. 

But even two-loop matching is mysterious. On general grounds, the classical string predicts the dimension 
of the operator to be \/AF(A/L^) with some function F{X/L'^) — L/VA + Any stringy nonpcrturbative 
contribution to the dimension should be of the form 

f{X/L')eM^^9WL'))- (1-3) 

Of course this term and others like it cannot be in general separated in the same exponential form on the gauge 
side, but just as an illustration let us momentarily take them seriously, in which case they can be perturbative in 
the gauge theory. In order to have a sensible contribution to perturbative gauge theory, we should require that 
the Taylor expansion of /(A/L^) and g{X/L^) about X/ = have integer and half integer powers respectively. 
If the lowest power of 5 is n < — ^ , then there would be an exponential suppression in L for all orders of A and 
H1.3|l would not contribute to any order in perturbation theory in the large L limit. However, if the first power is 
n = i, then H1.3|) can affect the perturbative expansion in the large L limit. However, H1.3|l should not contribute 
to the bare dimension of the operator, requiring that /(O) = 0. Hence, the term that can affect the perturbative 
gauge theory expansion at the lowest possible order would have the form ~ (A/L^)(l — aX/L + ^a^X^ /2L^ + ...). 
The lowest term in this expansion is a one-loop term, but it is suppressed by a factor of L from the classical 
string contribution so it can be ignored. The second term is a two loop term which scales like the classical 
contribution and so its existence would affect the two loop gauge contribution. The third term is a three-loop 
term which does not even have perturbative BMN scaling. What is puzzling is that the breakdown described 
here seems to take place at one loop higher. In fact it was shown in j42| that assuming that perturbative 
Yang-Mills is described by an Inozcmtsev chain ]3B], perturbative BMN scaling breaks down at four loops. 

We will demonstrate that indeed the Riemann-Hilbert problem on the string theory side immediately repro- 
duces one and two loops in gauge theory. We do not study the three-loop reduction of the sigma-modcl in this 
paper, but are planning ot return to this question in the near future. 

In what follows we will consider generic solutions to the class of so called finite-gap potentials. These 
solutions are directly related to algebraic geometry on finite genus Riemann surfaces. From the point of view of 
the Bethe equations, this corresponds to the Bethe roots condensing onto a finite number of "supports" in the 
complex plane. Such geometric objects have been previously used in SUSY gauge theories. For example, the 
Seiberg-Witten solution to J\f = 2 SYM is formulated in terms of the geometry of curves 0^1 which in turn can 
be nicely translated into the language of integrable systems ?505 in a class closely related to those considered 
below. In the Af — 2 case the trajectories of the integrable model correspond to BPS states while in the A/" = 4 
case they correspond to composite operators. This relationship deserves further investigation. 

We should mention that integrable spin chains arise in large- QCD [SH |52 1^ (the quasiclassical 
solutions to the Bethe anzatz equations were considered in this context in j5(il 157^ and in less supersymmetric 
cousins of Af — 4 SYM [HHl EHl EDI 1^ E21 • Potential relevance of the XXZ spin chain for gauge/string theory 
duality was discussed in 

In sect. 12 we review the Bethe ansatz and discuss the reduction of the Bethe equations in the scaling limit 
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to integral equations, similar to those found for matrix models. In sect.|31we construct the generic solution of 
this Bethe ansatz integral equation describing the one-loop A/" = 4 Yang-Mills theory and then generalize it to 
two loops. Sect. 21 is devoted to classical integrability: we show that the solutions of the Bethe equations in the 
scaling limit are in one-to-one correspondence with the finite-gap solutions of the weak-coupling sigma-model at 
one loop in the gauge coupling. Moreover, we then derive the classical counterpart of the Bethe equations in the 
string sigma-model and discuss the general underlying geometry of the sigma-model solutions. In particular, we 
find that the resolvent of the Bethe roots, which is the generator of the higher charges of the integrable system, 
is closely related to the quasi-momentum of the sigma-model and that the latter reduces to the former in the 
limit of zero coupling, i.e. at one loop. At two loops, the precise map between the sigma-model and the spin 
chain involves redefinition of the spectral parameter. In section |21 we apply the general formalism of sections 13 
and 0] to concrete examples. We obtain some new solutions as well as rederive known results to illustrate the 
general method. Included in these examples are general rational solutions where the Bethe roots condense onto 
one cut. We also consider a sigma-model solution which corresponds to a "pulsating string" and show that the 
quasi-momentum reduces to the resolvent of a particular solution of the £"0(6) chain in the scaling limit. In 
sect.iniwe present conclusions and speculations. Some technical issues are presented in Appendices. 



2 The Bethe Ansatz 

Here we will derive the Bethe-type equations for the diagonalisation of the one-loop dilaton operator. The 
generalisation to two loops will be presented at the end of the section 3. 

If we assign i?-charges (1, 0) and (0, 1) to the complex scalars €>i and $2, the eigenstates of the Heisenberg 
Hamiltonian with J down spins out of L total spins correspond to conformal operators with _R-charges ( Ji, J2) — 
(L- J,J): 

O = Tr + permutations) . (2.1) 

Each eigenstate of the dilatation operator (|1.2|l describes a collection of J interacting spin waves with rapidities 
ui, . . . ,uj. The spin waves attract and can form bound states, so the rapidities are in general complex. The 
eigenstates form representations of SU{2) generated by the total spin which commutes with the Hamiltonian. 
The highest weight in each SU{2) representation can be constructed with the help of the algebraic Bethe ansatz 
[IM\ and is completely characterized by a set of rapidities. These rapidities must be distinct complex numbers, 
are either real or form complex conjugate pairs, and satisfy a set of algebraic Bethe equations. 

The integrable structure of the Heisenberg spin chain is encoded in the transfer matrix (see, for example, 
pn] and references therein) defined as 



T{u) = Tr 



u+-<TL^cr\...\u+-(Ti 



(2.2) 



The Pauli matrices without a superscript cr act in a two-dimensional auxiliary space, and tracing in (|2.2I) 
over the auxiliary space leaves an operator in the Hilbert space of the spin chain. The transfer matrices at 
different values of the spectral parameter u commute and therefore generate an infinite set of conserved charges 
in involution ^. The Taylor expansion of the transfer matrix at u = i/2 generates mutually commuting local 

^To be more precise, for the chain of length L there axe exactly L independent charges, since the transfer matrix is a polynomial 
of order L. 
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charges: 

f{u) = (u+^] Uexp 



CO , 



(2-3) 

where U = e'^ is the shift operator that generates translations by one site. The first charge from the sum 
in the exponent of (|2.3|l is proportional to the one-loop dilatation operator H1.2(l : D = L + Q\ + O(A^). 
The Laurent expansion of p. 2(1 at u = oo (before taking the trace over the auxiliary matrices) produces the 
total spin and non-local Yangian charges (see the recent discussion of Yangian symmetry and its relation to 
Yang-Mills theory in 45 ). 

The Bethe ansatz diagonalizes the whole tranfer matrix, 

f{u)\u^...uj)^T(u)\u^...uj), (2.4) 
giving eigenvectors that depend on J rapidities and whose eigenvalues are 

* \ TT U — Uk — i I i \ -r-r U — Uk + i 



The second term does not affect the local charges (|2.3|l as it is multiplied by a huge power oi u — i/2. 

The Bethe equations are obtained from H2.5|l by simply noting that the function T{u), by construction, must 
be a polynomial of degree L, or for L — > oo an entire analytic function of u. Hence it has no poles at any u — Uj, 
and this immediately leads to the J conditions 

^V= n (2.6) 



which are in fact the Bethe equations. 

These Bethe equations also represent periodicity conditions for the wave function and can be written in the 
logarithmic form as 

/ Ui + il2\ v-^ — Uh + i 
L log ^—^ = J2 — - ^'^^^J - (2-7) 

The mode numbers rij are arbitrary integers. The cyclicity of the trace in the operators of A/" = 4 supersymmetric 
Yang-Mills requires that 

4P -TT Ui + i/2 , ^ 

^^"-n^=i' (2-8) 

which is the quantization condition on the total momentum: ^ G The energy of the solution gives the 
anomalous dimension 7 of the corresponding operator in Yang-Mills theory |27| , 

j=i J 
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In the scaling limit of large operators where L ^ oo, all Bethe roots are of order uj ^ L. By rescaling Uj = Lxj 
we can then write eq. H2.7() in the form 

^^I ^ ^^-2™,- (2.10) 

Let us consider the situation with a finite number of different mode numbers rij , and assume that the number 
of Bethe roots with the same mode number is of order L. The roots with mode number are centered roughly 
around the point x = 1/(27™^), with a typical separation between adjacent roots of order Ax ^ 1/L. Hence, 
roots with the same mode number form a continuous contour in the complex plane in the scaling limit L oo^. 

Equations (|2.1(J|) are very similar to the saddle point equations for matrix models (apart from the last term, 
to be discussed below), where the Bethe roots are the analogs of the matrix eigenvalues. Just like for matrix 
eigenvalues, one can characterize the distribution of Bethe roots by the density 



P{x) = j^S{x-Xj) 



L 



(2.11) 



or by the resolvent 



L f-^ x-Xj Jc X- E, 



The density is non-zero on a set of contours C in the complex plane which in general consists of several non- 
overlapping cuts, C = Ci U . . . U Cki where the Z**^ cut C; — {ai,bi) represents roots with mode number m. 
Though we will assume that the number of cuts is finite, this requirement can certainly be relaxed. We will 
return to this point in the discussion of classical integrability of the string sigma-model. As follows from its 
definition, the resolvent is an analytic function on the complex plane with cuts. Asymptotically it has the form 

G{x)^ - + ..., (x^oo), (2.13) 

X 

which reflects the normalization condition for the density 

d^p(0=a, (2.14) 

c 

where a = J/L is the filling fraction, the fraction of down spins in the Bethe state. Strictly speaking, Bethe 
states only make sense for a < ^ , but it was observed in a number of examples |32[ I33L that solutions of 
the Bethe equations can be analytically continued past a = ^. We will clarify the meaning of this analytic 
continuation in sections l51and l5T^ The analytic continuation is always possible and the analytically continued 
solution corresponds to some Bethe state with a replaced by 1 — a. Since replacing a with 1 — a is equivalent 
to interchanging $i and $2 in the operator Tr($^^$2^ + ••■)' = L — J, J2 ~ J, the resulting state has 
essentially the same charges. We thus allow for < a < 1 and regard charges of the states with a and 1 — a as 
equal. 

•^The limit we consider here is different from the standard thermodynamic limit. We are keeping mode numbers finite, while 
in the thermodynamic limit, momenta pi = 2-Kni/L are fixed as L goes to infinity. Our Bethe equations are reminiscent of those 
for the so called Gaudin model. 
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The scaling limit of the Bcthc equations H2.10|l can be rewritten as an integral equation for the density, 



G{x + iO) + G{x - iO) ^ 2J ^ 1 + 27rn;, x e Ci (2.15) 

J c x-^ X 

The momentum condition (|2.8() after rescaling turns into 

1 ^ 1 

= 27rm, meZ (2.16) 



j=l J 

or, taking into account the Bethe equations (|2.1()(l and H2.15|l . 



L / ^(^)^^ = V / Pi^)d.^' ^ 27r Vn; / p(x)dx = 2TTm, ni,me'Z- (2-17) 

27rz/c a; ^1 2; ^ 7c, 

For the anomalous dimension (|2.9|l we obtain in the scaling limit 

Rescaling the spectral parameter u = Lx + i/2 in the eigenvalue of the transfer matrix and dropping the second 
irrelevant term in (|2.5() we find in the limit L —^ 00: 

T{Lx + i/2) _ i ^ ST ^ \ 

[Lx + il2t '"''^^yi.^x^ x,j ^ ' 

Comparing with H2.3|l . we see that the resolvent (|2.12l) is the generating function of local charges in the scaling 
limit 021 1^ . Obviously, the first two coefficients in its Taylor expansion are the total momentum and the 
energy. ^ 

Our derivation of the integral equation for the density of Bethe roots rested on the assumption that the 
right-hand side of the Bethe equations can be expanded in powers of l/{uj — Uk)- This is not always true. The 
assumption breaks down if the distance between nearby roots Uj — Uk approaches i. Some solutions of the Bethe 
equations are known to have roots separated exactly by i in the thermodynamic limit. The roots form arrays 
with equidistant separation: Uj = uq + i j. Such arrays are usually referred to as "strings" |31| and survive 
the thermodynamic limit due to the chain cancellation of dangerous terms in the sum entering eq. H2.7|l . It is 
certainly possible to have large strings composed of a macroscopic number of roots. The density will then have 
constant pieces with a flat distribution, p{x) = 1. Such condensates of Bethe roots necessarily have tails with 
smaller density. Condensates thus can be viewed as extra logarithmic cuts in the complex plane that connect 
pairs of the original contours C; (the situation is identical to that arising in multicolor 2D Yang-Mills theory 
on the sphere ^64^). The resolvent experiences a jump by 2t: across any of these extra cuts and is no longer a 

^There are only L independent charges, so only the first L coefficients of the T(u) Taylor expansion make sense as conserved 
charges. Unless we are interested in non-local charges, we can forget about the second term in 12. 5i . The generating function of the 
local charges is a pure phase: exp (jHermitean operator) and not a cosine of it. In sect. |3 we get a cosine of a Hermitean operator 
because T{u) is expanded at the "wrong" point u = (instead of m = i/2). In effect, the argument of the cosine is shifted by 
compared to the argument of the exponential. 
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single- valued function. The integral equation (|2.15|) stays the same, but we have extra conditions on its solution 
as it crosses these logarithmic cuts. 

The solution of the Bethe equations with a condensate along the imaginary axis is one of the examples in 
which the agreement between the anomalous dimension and the energy of a string soliton was first observed 
j32j . In that example, however, the roots were separated not by z, but by a half-integer interval i/2, and the 
condensate had the density p{x) = 2. We do not see any reason why strings of an arbitrary integral density 
cannot appear: for the distribution of roots Uj — uq + i where m is any positive integer, the same cancellation 
mechanism is at work |ti5| . Though these denser strings apparently have never been discussed in the literature 
on the Bethe Ansatz, their existence is quite natural, at least in the long chain limit considered here. Another 
argument in favor of condensates of an arbitrary density is that solutions p{x) — 2m are necessary to match 
the string theory predictions, where the integer m plays the role of the winding number. 

To summarize, Bethe roots in the scaling limit are continuously distributed on a collection of contours C; 
in the complex plane and may also form condensates, which connect pairs of contours and on which the density 
takes constant integer values. The density of roots p{x)dx is real, positive and is normalized to the filling 
fraction a. It must also be symmetric under complex conjugation of x. The density of Bethe roots satisfies an 
integral equation H2.15|l which holds on each of the C;'s. All possible solutions to the Bethe equations in the 
scaling limit are thus parameterized by specifying mode numbers n/ associated with contours C; and integer 
densitities of condensates to/. Once these numbers are specified we can solve the integral equation, at least in 
principle, and find the density. We will see, however, that there is a symmetry that transforms different sets of 
integers (n;, to;) into each other so that solutions with different sets of to's and n's lead to the same resolvent. 

3 General solution 

In this section we first construct the general solution of equation 1)2. 15|) and then consider corrections to this 
solution arising from the two-loop contributions in SYM theory. 

3.1 General one-loop solution 

Following the procedure proposed in we introduce 

p{x) = G{x) - ^, (3.1) 
2x 

which we shall call the quasi-momentum for reasons explained in the next section, and rewrite 1)2.15(1 in the 
form of a Riemann-Hilbert problem defining the analytic function p{x) from its discontinuities 

p{x + iO) + p{x - iO) = 27rni, x ^ Ci (3.2) 

on every cut C^. Equation ((3.2(1 together with relation 1(3.1(1 are very similar to the corresponding equations 
arising in the quasiclassical computation of matrix integrals. The main difference is the presence of the mode 
numbers n^, which are different on the different eigenvalue supports. Therefore, unless there is only one cut, 
the resolvent G{x) and the quasimomentum p{x) are not algebraic functions, i.e. in contrast to a very similar 
equation found for the one matrix model (see e.g. [HSj) and even the more complicated situation of the two- 
matrix model '67' , the resolvent and quasimomentum can not be determined by solving a polynomial equation ^ . 

^This is a crucial diflference between the quasiclassical solutions to the Bethe equations we consider below, as compared to those 
of I5KI . which correspond to the limit J ^ oo aX, finite L. In the last case the resolvent and quasimomentum are logarithms of 
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It follows nevertheless from ()r>.2|l that p{x) is a double- valued function of x. It must be single-valued on the 
physical sheet of its Riemann surface, which is a complex plane with cuts Cj. The only singularities of p{x) are 
two poles at a: = on the two sheets of the Riemann surface. The pole on the physical sheet should cancel in 
G{x), which by definition should be regular at zero. In other words, p{x) = J dp is an Abelian integral for the 
meromorphic differential dp, which has two double-poles at a; = and integer periods on a hyperelliptic curve 
S. Let 

2K 

S: y^ = R2K{x) = x^'' + nx^''-^ + ... + r2K = Y[{x-Xj) (3.3) 

where the product is taken over the end-points of all K cuts C; e C. We should also impose the reality 
constraint on the whole solution such that the branch-points of every cut have complex conjugated coordinates, 
= y^2j+i- The other possiblity is that some x^'s are real, but it is never realized in the present case, although 
solutions with real cuts do exist for spin s = chains and are relevant for SYM duals of strings with 
nonzero angular momentum in AdS^ In another context the quasiclassical solutions to Bethe equations 

with the cuts along the real axis were discussed in |56|. 

Any meromorphic differential on (|3.3|) has general structure df{x) + Sisl^^ where f{x) and g{x) are two 
rational functions. From (|3.2I) one can expect for that differential dp is just dp — g{x) and we can rewrite 
H3.2|l in terms of integrals over the B-cycles on E, 

27r(7ii ~ rij) = p{xi + iO) — p{xj — iO) + p{;Xi — iO) — p{xj + iO) = 

rx,+iO rx.-iO . (3.4) 

= / dp + dp = (p dp. 

Jxj-iO Jxj+iO J Bij 

Here we used the hyperelliptic involution y ^ —y of the Riemann surface S to move the contour of integration 
from the physical sheet to the second sheet. The involution interchanges the two sheets and flips the sign of 
y, then the differential dp changes sign too and hence we should simultaneously reverse the orientation of the 
contour, this is illustrated in fig. ^ 

algebraic functions and there is no freedom in choosing discontinuties in lij.2l . while below, together with the filling fractions of 
eigenvalues, these will be the parameters of our general solution. 
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Figure 2: Hyperelliptic Riemann surface S, defined by equation (|3.3() : zig-zag line represents the condensate 
cut. There is a double pole of differential dG at the image of the origin a: = on the second sheet, while the 
differential dp has double poles at a; = on both sheets. The sheet where dG has no singularities is called 
physical. 

Equation (|3.4|) imposes an integrality condition on the B-periods of dp, 

dp = 2Tr{ni - uk) i = l,...,K -1, (3.5) 

for the canonical choice of the B-cycles on the hyperelliptic surface, Bi = Bm (see fig. EJ- In the absence of 
condensates the single- valuedness of p{x) on a physical sheet requires in addition that A-periods are zero, 

j dp = i = l,...,K-l, (3.6) 

for the canonical Ai cycles on the hyperelliptic surface (fig.EJ. 

We need an extra condition to satisfy all K equations in (|3.2() . To find it, one can rewrite the K-ih equation 
in as 

p(oo+) -p(oo_) = / dp = 2TmK, (3.7) 

J CxD_ 

where by oo+ and oo_ we denoted two points with a; = oo on the first and second sheets of S respectively. The 
contour of integration comes from the infinity on the physical sheet, passes through the K-th cut Ck onto the 
second sheet, and then goes to the infinity on this sheet (fig. OJ. The points oo+ and oo_ are marked by the 
fact that the differential xdp, which determines the number of Bethe roots, has first order poles there (with 
residues of opposite sign), though dp itself is non-singular. Note that since the choice of the cut K is arbitrary 
we could choose the set of equations similar to (|3.7|l on all K cuts to replace the equations H3.5|l . 
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Figure 3: The contour of integration in H3.7(l . The marked points on two sheets correspond to infinities 00+ 
and oo_ where x = 00. 



Since dp has second-order poles at a; = on both sheets of the Riemann surface, it is of the general form 

A'-l 



, dx \ ^ 

dp = — > 

The pole at x = on the physical sheet must be cancelled in G{x), which requires dp = 
the first two coefficients to be^ 



akx 



a-i 



/r2K 



ao 



r2K-i 



(3.8) 

0(1) and fixes 
(3.9) 



The system of equations (|3.6|l uniquely determines the rest of the K — 1 coefficients and thus completely 
defines the differential (|3.8(l on a given Riemann surface 1)3.3(1 . The rest of the conditions are conditions on the 
moduli of the Riemann surface itself. 

All together the curve 1(3. 3|l is parameterized by 2K end-points of the cuts . The integrality of B-periods 
expressed by the eg s . ((3.511 and the extra condition 1(3.7(1 constitute a system of K equations on the coefficients 
of 1(3.3(1 and thus reduce the total number of free parameters to K. These free parameters exactly correspond 
to the filling fractions of the Bethe roots on the K cuts Cj , 



S3 ^ 77- f P{x)dx = [ p{. 



x)dx, J ' = 1, . . . , if — 1 



Sk 



1 
27ri 



p{x)dx 



p{x)dx 



a 



Ck 



i=i 



(3.10) 



^We assume that the branch of the square root is chosen such that y = x 



at infinity on the physical sheet. The value of 



y at zero may coincide with the algebraic square root of r2K or may differ from it by a sign. The square root ^r2K in this formula 
should be understood in the analytic sense as yifS). 



11 




The residue at infinity determines the total number of Bethe roots ^: 

i - a = --^reSooS^dp (3.11) 
1 Am 

The normalization condition immediately allows us to set aK-\ — \ — a. The total momentum H2.17|l can be 
rewritten in terms of the filling fractions Sj in H3.1()(l as 

P = ^reS:j=o ^i^ll^ = 27r V n.S. = 27rm. (3.12) 

The anomalous dimension H2.18|l can be now expressed directly through the parameters of the hyperelliptic 
curve (|3.3|l and of the quasimomentmn 13.8|l or resolvent G{x). From the expansion of p(x) at x = 0, 

, , 1 „ S-K^Lj 2n 
P[X) = ZTTm X + Uix ), 

2x A 

we have 



o 1 9 I (3-13) 

Stt^L V 4r2K ISr-a/f \fr2K ) 

which is a general solution for the one-loop anomalous dimension of "long" operators (|1.1|) . It is expressed 
explicitly through the coefficients of the equation (|3.3|) which defines the complex curve, and by which is 
determined through the solution of the linear system H3.6|l . and is a transcendental function of the coefficients 
of (|3.3() . The coefficients of (|3.3|) are themselves defined implicitly as functions of the mode numbers rij, the 
total momentum P = l-nm^ and the filling fractions Sj. 

Very little changes if we allow for condensates. The differential of the quasi-momentum is still an Abelian 
differential of the second kind, the equations H3.2|l are still solved by imposing integrality of B-periods on d-p 
and an extra condition l|3.7|l . The only difference is in the A-periods. If an A-cycle Ai crosses the condensate 

''^Note, that 13.111 directly corresponds to the asymptotics t'2.131 
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cut, p jumps by 2wini, where rrii is an integer. Hence, A-periods are no longer zero but are also integral (fig-EJ. 
We thus get a very symmetric set of conditions on the periods of dp, 



dp = 27r(ni — nx), 

(3.14) 

dp = 2'Kmi 

A. 

The counting of parameters does not change. There are 2K — 1 equations (A and B cycle conditions and (|3.7I) ) 
on the total of 2>K — 1 parameters {2K moduli of the hyperelliptic curve S and K — 1 free coefficients in the 
definition of dp. eq. H3.8|) '). The K parameter freedom reflects the possibility to redistribute Bethe roots over 
the K cuts with arbitrary filling fractions (|3.1U|I . The overall normalization H3.11|l fixes the total number of the 
roots and leaves K — 1 independent parameters. 

The approach of this section allows one to calculate the resolvent and consequently all the conserved charges 
for a given Bethe state in the scaling limit. It may seem that Bethe states are in one-to-one correspondence with 
the sets ofK filling fractions and 2A'— 1 integers {rii, rrii], but this is not quite true. Not all sets of parameters are 
allowed, since some of them give rise to complex curves that do not satisfy the reality condition, which therefore 
is a constraint on allowed values of integers and m^. Also, some different sets of n's and m's may correspond 
to the same Bethe state. Suppose that we solved for the resolvent and found all its branch points Xj as functions 
of the parameters. There are still many ways to connect these branch points by cuts C/. In principle, the cuts 
should be defined in such a way that the density p{x)dx for p{x) = DiscG(a::) = {G{x + iO) — G{x — iO)) 
is real and positive definite on the cuts. Near each branch point dx p{x) ^ dx {x — Xj Y^^. There are three 
lines on which dx {x — real and positive, and therefore three cuts consistent with positivity can meet 

at each of the branch points. This leads to a discrete ambiguity in cutting the complex plane. This argument 
is standard in the discussion of multi-cut solutions of matrix models (see, e.g. j6()l I68p . Any rearrangement of 
cuts or permutation of the branch points induces a linear transformation of the A- and B- cycles and therefore, 
according to H3.14(l . of the integers Ui and m^. The branch of the square root on the physical sheet may, however, 
change sign at infinity under such rearrangement. The change in sign of the residue of xdp at infinity, according 
to (j3.11|l . is equivalent to interchanging a with 1 — a. If a > ^ from the outset, by rearranging cuts we can 
always get a new a in the physical domain < a < i. Therefore, an analytic continuation of a solution of the 
Bethe equations beyond a = ^ indeed corresponds to a physical Bethe state with a < ^. The energy and all 
local charges are modular invariant and do not change under the rearrangements of cuts. Since local charges 
completely characterize a Bethe state, we can relax even the reality condition for the density, and connect the 
branch points by cuts in an arbitrary way. Then all Sp{K — 1,Z) transformations of cycles are allowed. In 
particular, it is always possible to set all A-periods of dp to zero by an appropriate Sp{K — 1, 2) transformation. 
In other words, it is always possible to get rid of the condensates. This remark will be very important in the next 
section, when we will compare general solutions of the Bethe ansatz with classical solutions of the sigma-model. 
This discussion may seem somewhat abstract, but we hope to clarify it with explicit calculations for two-cut 
solutions in sectional 

A similar ambiguity affects the total momentum H3.12|l 

P = -G(0) = / ^ dG. (3.15) 
Jo 

The contour of integration connects zero and infinity on the physical sheet and is otherwise arbitrary. Adding 
an A-cycle Ai to it shifts the momentum by 27rmi, but the momentum is defined only up to an integer multiple 
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of 27r, and only e'^ makes sense. The case when the condensate passes through zero [32] requires special care, 
because in that case the integral in (|3.15() is ambiguous and should be defined symmetrically such that the 
corresponding A-cycle contributes half of the period into the momentum. With this prescription the condensate 
of density m (Bethe roots stretched along the imaginary axis with the spacing im) has the momentum mn. 
Indeed, careful inspection of the general formula 12.8|l shows that an evenly spaced condensate has e*^ = 1, 
while an oddly spaced distribution of roots gives e*^ = —1. 

The Bethe equations H2.7II have an obvious symmetry, Uj —> ~Uj, in which case the numerator replaces the 
denominator and vice versa. Therefore, sometimes it is rather natural to consider a class of solutions, when 
magnetic numbers satisfy n(—Uj) = —n{uj), Then the curve (|3.3fl has a symmetry x —x, i.e. 

K 

2/2 = [](x2-x2) (3.16) 

The counting here is slightly different. We start from K independent parameters (instead of 2K), impose K/2 
constraints (|3.7() . since K should be even in this case, and the rest of parameters are eaten by the residue of 
the differential xdp (|3.11|l and half of the independent "fractions" H3.10|l . A particular example of a symmetric 
solution will be considered in section [^31 



3.2 Two-loop corrections 

By this approach we can find also the two-loop corrections to our general one-loop solution, thus diagonalizing 
the operator (|1.2|) in the classical limit, using the two- loop perturbative integrability hypothesis proposed in 
and confirmed in [691 EU] . A way to take higher loop corrections into account is to consider the Inozemtsev 
generalization of the spin-chain in the hyperbolic limit j48| . as was recently proposed by Serban and Staudacher 

m 

The Inozemtsev chain has long range interactions between the sites, refiecting the fact that g-loop corrections 
to the dilatation operator involve spin-exchange separated by up to g sites. The hyperbolic limit has one free 
parameter that describes the fall-off of the interaction and Serban and Staudacher |42| mapped this parameter 
to X/L^. They then showed that not only did the Inozemtsev chain reproduce the two and three-loop dilatation 
operators in [32112^1 EDI j but it also matches the two-loop predictions for the string solitons in [201122]. However, 
the three-loop predictions in [32] differ with those in |2l)[ I22|. and furthermore, there is an explicit violation 
of BMN scaling starting at four loops [32] ■ Despite the fact that the equivalence of the full perturbative SYM 
to the Inozemtsev chain is still only a hypothesis, the two-loop corrections are robust, and so we can use it to 
study the anomalous dimensions of general operators at the two-loop level and compare these with predictions 
of quasiclassical string theory. 

The two-loop corrected Bethe equations look as follows 02] 



fc=l(fc#j) 

where the momentum pj — pj(uj) is defined by solving the equation 

A 



n ^^i^^, (3.17) 



u = uo{p) + -^smp + Oi\^) (3.18) 
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with the function uq(p) given by 



uq + i/2 
uq — i/2' 



Up to two- loop order, eq. H,S.17|I can be rewritten as 

Uj + i/2 



87r2 u2 + l/4 



n 



Uj — Uk + i 



"^3 V2 8;r2 u^.+l/i 

\ 3 ' , 

This means that in the classical limit u = xL, L ^ oo, eq. H2.15|l is corrected as follows 

2-f x&Ci 



(3.19) 



(3.20) 



(3.21) 



with the same normalization of the density (|2.14|l as in the one-loop case. Note however, that being perturbative 
from the gauge theory side, the second term in the l.h.s. of H3.21|l is "singular" perturbation of the integrable 
system, since it adds a singularity of higher degree. For example, analogous perturbation in the algebraic matrix 
model case jfifil 16 7| would change the genus of the curve. Here, it will causes some problems, when comparing 
with the string theory side. 

The two-loop generalization of the total momentum quantization condition H2.8() can be rewritten as 



n 



-1 \ Uj — i/2 



= 1 



or, in the classical limit, 



27rm 



dx p(x) 



1 A 1 

-I- 



x Stt^L^ x^ 
The anomalus dimension is corrected as follows 

J 



dx 



,G{x) 



A 



J 

E 



A 



1 



• 2 

— ^ sm-^ 
27r2 2 



3A2 



) A 



2m 



^sin^^ 
87r4 2 



A 1 



87r2 1*2 + 1/4 1287r4 (u^ + 1/4)2 i287r4 {u^ + 1 /4)3 



0(A3) 



In the classical limit L ^ oo this two- loop formula corrects ea. H2.18() in the following way 



^^Ljdxpix) (^^ + 



3A2 



1 



1287r4L4 a;4 
3A2 



^4 



A 1 



87r2L a;2 1287r4i3 ^ 



^ + 0(^)+0(A»)) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



We kept here only the terms relevant in the BMN limit of fixed T 
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We can now solve eq. (13.2111 in the same way as (|2.15|l . We introduce 

p{x) = G(x) ^^--^ (3.26) 
and find for the differential dp on the Rieniann surface (|3.3|) the following modification to H3.8|l 

dp=—y akx''-\ (3.27) 

^ fc=-3 

with the condition that dp = + ^ + 0(1) for x ^ 0. It is clear that this gives two extra conditions to fix 
the two extra coefficients a_3 and a_2- 

Below we will compare these results with the two-loop correction following from the sigma-model results, 
where the term ^ will appear from the two-pole structure \ {^-^^ + — ^ ~^ ^ ~^ ^ {^) integral 

equation similar to H3.21|l . We will see that although the structures are very similar, the details are different. In 
contrast to one loop, where the equations are easily reproduced from those of the string sigma model, the two- 
loop equations of this section follow from the integrable equations of string theory only after a reparameterization 
of the spectral parameter. This basically reflects the fact that the choice of the spectral parameter in (|3.18(l 
and H3.21|l is somewhat arbitrary. Nevertheless, they lead to the same results for the anomalous dimension up 
to the two-loop level. In sectional we will present a wide class of examples where the agreement is explicitly 
checked. 



4 Classical integrability 

We now turn to integrability on the string side of the AdS/CFT duality. Classical integrability of the sigma- 
model in the ^^5*5 x background ^503| means, among other things, that the classical equations of motion 
of a string in AdS^ x can be solved by algebro-geometric methods (HHl IZOl [31 IZ2] or, for certain particular 
cases, by the inverse scattering transformation |73j . In this section, we will examine finite-gap periodic solutions 
of the sigma-model and establish direct links between them and the Bethe ansatz. A pedagogic introduction 
into finite-gap integration methods can be found, for example, in the monographs j71[l72| . and the comparison 
of quantum inverse scattering with classical integrability, applied mostly to the nonlinear Schrodinger and KdV 
equations, can be found in 36 , 37^ . 

As in the discussion of the SYM operators we will focus on an SU (2) reduction of the full sigma-model to 
the subsector of string moving on x Rf. Many explicit solutions are known in this case |26|. The string 
action in the conformal gauge is 



So 







/ da 


[ dr 


Jo J 





(4.1) 



where Xq is the global AdS time and Xi, i = 1, . . . ,4 are Cartesian coordinated on embedded in R^: 

X,X, = 1. (4.2) 

All other world-sheet coordinates (one radial and three angular coordinates on AdS^ and two extra angles on 
S^) are set to constant values. Such an ansatz describes a string localized in the centre of AdS and moving in 
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an subspace of 5^. The effective string tension is related to the 't Hooft couphng according to AdS/CFT 
|7]. The equations of motion that follow from 14.1|l should be supplemented with the Virasoro constraints: 

{d±X,)' = {d±Xo'f - k', (4.3) 

where a± — ^(r ± a), d± — dr i da, and we always consider a gauge Xq — nr. 

It was shown in [^Hl that the equations of motion for the string whose centre of mass moves along a big 
circle of with a large angular momentum have a well-defined weak-coupling limit (A/L^) <C 1, where now 
■\/A is the effective string tension and L is the angular momentum, and reduce in this limit to^ 

smedrO + 9^(sin2 ed„ip) = 

(4.4) 

die ^ sinedr'f - sine cose{d„Lpf ==0, 

where ip and 0, in addition to another angle that parameterizes the motion of the centre of mass, are coordinates 
on . This limit can be regarded as an approximation of an almost light-like string motion |74| In terms of a 
unit vector, S — (sin0 cos (p, sin sin (p, cos 0) the equations (|4.4(l take the form 

= S X S,,. (4.5) 

This is the equation of the classical Heisenberg model. It can be formally obtained from operator equations of 
motion in the spin chain by replacing quantum spins cr; with C-number unit vectors S and simultaneously taking 
the continuum limit: <t/(t) — > S{2ttI/ L,t) — S((t, r). In |35) Kruczenski presented a path-integral derivation of 
this fact, and identified this classical approximation with the limit of the large size of the chain in which only 
slowly varying spin configurations are kept. 

The classical Heisenberg equation is completely integrable [TSJ |7B1 [73 [73] and is a particular case of a more 
general Landau-Lifshitz equation 

St = S X S^^ + S X JS (4.6) 

where J is an anisotropy matrix which is proportional to the unit matrix in our case, so that second term in 
the r.h.s. vanishes and (|4.6|) turns into 1)4. 5|) ^. We examine the classical Heisenberg equations in the next 
section and establish an equivalence between their periodic solutions and the general scaling solution of the 
Bethe ansatz which was obtained in section [2| 



4.1 Finite-gap solutions of the Heisenberg magnet 

In order to relate integrable structures of the classical and quantum Heisenberg models, let us perform a 
continuous limit directly to the transfer matrix l|2.2|l . If we replace cr; in (|2.2(l by S{2Trl/L) and also rescale the 
spectral parameter u — Lx, then the transfer matrix becomes 



T{Lx) 



n (l + ^ S{27Tl/L) <t) 



**Our normalization of the world-sheet coordinates t and a is different from that of 1351 . We take < cr < 27r and rescale t by 
the total spin L, consistently with the normalization S^ = 1. 

^Note that the dynamical Neumann systems considered in 13"- 26 can be considered as a particular reduction of equation 14. 5i . 
and we do not necessarily need to introduce the anisotropy 14. 61 . contrary to the proposal of 1631 . 
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This expression has a weh-defined continuum Umit, as an anti-path-ordered exponential 



rifl + ^^(27rVi)®^) nix)^Pexp(^<f^S{a)-cT) (4.7) 



Tiflix) = hm ' = 2cosp(a;) (4.8) 



which is just the monodromy matrix of the classical Heisenberg model, and we can introduce the classical 
quasimomentum p(x) by 

T(Lx) 

L^oo {Lx) 

The classical equations of the Heisenberg magnet are equivalent to the consistency condition [L, M\ = of an 
auxiliary linear problem 

£»= (a, -^)»(i;t,,t) = 



where the Lax matrices are 

[/ ^ cr • S 
F = cr • (S^ X S) 



(4.9) 



(4.10) 



The monodromy matrix is a parallel transporter of the flat connection (L, M ) around the circle: ^(x) — 
^'(a;; tq, 1), where ^'(x; r, a) is the solution of 1)4. 9|) with the initial condition '^{x\ tq, 0) = 1. Since the connection 
is flat the trace of the monodromy matrix does not depend on tq and therefore generates an infinite set of 
integrals of motion. The classical monodromy matrix is unimodular and unitary when the spectral paramater 
is real. Its eigenvalues e ^'p(^) then determine the quasi- momentum p{x) , and the Taylor expansion of the 
quasi-momentum generates an infinite set of conserved charges, which are local functionals of S(r, a) [73]. 

From section 121 we know that taking the scaling limit of the quantum transfer matrix in H2.5|) leads to 

-^ = 2cosfG(x)-i 
{Lx)^ \ ^ ^ 2x 



2 cosp(x) ^ = 2 cos ( G{x) - - ) , (4.11) 



with the resolvent defined in (|2.12l) . There is a l/(2x) correction term in (|4.11l) because, unlike in H2.3|l . we did 
not shift the spectral parameter by i/2. We conclude that the resolvent constructed from the density of Bethe 
roots and the quasi-momentum of the classical problem are related exactly as in (|3.1|l . The local conserved 
charges in the spin chain and in the classical Heisenberg model are thus equivalent. This observation elucidates 
the meaning of the transformations used in j46| to relate higher conserved charges in the sigma-model to the 
conserved charges of the spin chain. The approach of 021 is based on Backlund charges, which differ from 
charges defined by the expansion of quasi-momentum. One has to do a linear transformation on an infinite set 
of charges to account for this difference. 

The transformation of the quantum spin operators into the c-number unit vectors is hard to motivate 
rigorously in the present context. There is a more direct way to compare the Bethe ansatz of the quantum 
spin chain and the auxiliary linear problem of the classical model, which allows one to derive H3.1|l directly, 
by inspecting the analytic structure of the quasi-momentum. Moreover, this method can be generalized to the 
sigma-model, where the quantum counterpart (all-loop Bethe anzatz) is not yet known. 
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Figure 5: The trace of the monodromy matrix as a function of the spectral parameter along the curve on which 
cosp(a;) is real. The crosses correspond to degenerate periodic solutions, which represent forbidden zones shrunk 
to points. 



It is always possible to reconstruct the solution of the Heisenberg equation from the solution of the auxiliary 
linear problem. Indeed, it follows from (|4.9() . and H4.10() that 

S = -2TTix Tr {a- 1'^^'"^) . (4.12) 

A slightly less trivial statement is that ^'(a;;r, cr) is essentially determined by the analytic properties of the 
quasi- momentum p{x). Let us now establish these properties by inspecting the auxiliary problem 1)4.9(1 . 

The linear problem (|4.9() has a form of a one-dimensional Dirac equation 

were V' is a two-component spinor. Two linearly independent solutions of this equation can be chosen quasi- 
periodic. Indeed, if the initial conditions are eigenvalues of the monodromy matrix 

Vt{x)i}±{x-Q) = e±*f(^)V±(a;;0), (4.14) 

the solution 4'±{^'i^) = '^{x',<j)ip±{x]Q) will again satisfy 4'±{x]a + 27r) = e ^*P^^^?/'±(a;; cr) because 4'(x;cr -f- 
27r) = '^{x\(j)n[x). 

However, the quasi-momentum is not neccesarily real, being in general, a complex-valued function of the 
spectral parameter. The condition that its imaginary part vanish is a single equation on two variables (the real 
and imaginary parts of the complex spectral parameter x). Therefore, the quasi-momentum is real on a set 
of disjoint one-dimensional supports (some segments in the complex plane of x). These segments correspond 
to allowed zones of the spectral problem 1(4.13(1 . In a more familiar setup of the Schrodinger equation with a 
periodic potential, the allowed energy zones lie on the real axis, but the operator ((4.13(1 is not Hermitean. The 
allowed zones are singled out by the condition that Im (Tr fi) = and |Tr f7| < 2. One can also define forbidden 
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zones as loci where Im (Trfi) = and |Trf2| > 2, see fig. El In general the number of allowed and forbidden 
zones is infinite. We will only discuss finite-gap, or algebro-geometric solutions for which this number is finite. 
They are governed by the geometry of a complex curve of finite genus H3.3|l and directly correspond to solutions 
of the Bethe equations with a finite number of cuts. In other words, it is more convenient to work with the 
smooth non-degenerate curve ]3.'6\i than with the curve, defined by the Baxter equation on quasimomentum 
det{n - e^P) = or 

e'P + e-'P = Ti-nix) (4.15) 
which is generally a degenerate curve of an infinite genus^". 

At zone boundaries the monodromy matrix degenerates into the Jordan cell and has only one eigenvector 
with an eigenvalue 1 or —1. The quasi- momentum becomes an integer multiple of tt and the two independent 
solutions of (|4.13|) become degenerate. This degenerate solution of H4.13|l is either periodic or anti-periodic. 
Accidental periodic and anti-periodic solutions may appear inside some of the allowed zones (see fig.lSJ. They 
are always present for finite-gap potentials, and this is related to the fact that p{x) is not an algebraic function 
on (|3.3|l . The linear problem (|4.13() still has two independent (anti)-periodic solutions at these special points. 
Hence, is either -1-1 or —1, and these points do not correspond to zone boundaries. They rather represent 
forbidden zones shrunk to a point and of course are not generic, as much as finite-gap solutions are not generic 
in the space of all solutions to the Hcisenberg magnet equations. 

The zone boundaries are singular points for ■iP±{x;(t) as functions of x, because ip± collapses into one 
degenerate solution there. It is thus natural to regard and tp- as two branches of a single analytic function 
•0 on a double cover of the complex x plane. The zone boundaries are branch points at which the two branches 
collide. Apart from an essential singularity at a; = 0, V' is a holomorphic function on a hyperelliptic curve, two 
sheets of which are glued together along the forbidden zones. This is how the Riemann surface S, defined by 
1)3. 3|l . arises in the auxiliary problem of the classical Heisenberg model. The two eigenvalues of the monodromy 
matrix, e^'^*^^-* are also branches of a single meromorphic function on S. In fact, the trace of the monodromy 
matrix is an entire function of x. There is no reason for it to be singular anywhere, except at a; = where the 
potential in H4.13|l has a singularity. But solving H4.15|l for p we will encounter square root singularities when 
the discriminant of the quadratic equation H4.15|) . with respect to e'^, becomes zero. This is another way to 
see why the quasi-momentum has branch points. 

The matrix solution of H4.9|l . ^'(x;r, fi), is a meromorphic function on S as well. This latter function has 
the following asymptotic behavior: 

4'(x;T,cr) expfiffaf-^-K— 5-^)1 •$(a:;T,cr)-C (4.16) 



V47ra; An^x' 



where C is just a constant matrix (independent of x, r and a), and 

<^{x; T, a) = (j)o{T, <J) + ^ <Pk{T, a)x^ ^4 ^^y-j 

fc>o 

is regular at x = 0. An analytic matrix function 14.16|l on the hyperelliptic Riemann surface 1)3.3(1 with such 
asymptotics and g + \ extra poles, which can be chosen independent of r and ct, is unique, they can be expressed 
in terms of the Riemann theta-functions j78| by the standard technique of finite-gap integration. 

^"Note, that here is an essential difference with the quasiclassical chain with a finite number L of large spins. In the latter case 
14.151 basically coincides with 13. ill and the maximal number of cuts is bounded by the length L of the chain. We restrict ourselves 
to this wide class of finite-gap or algebro-geometric solutions, though the equivalence between the Bethe ansatz and the classical 
solutions certainly holds for infinite-zone solutions too. 
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Let us return to the analytic properties of the quasi-momentum. Consider the values of the quasi-momentum 
p{x + «0) and p{x — «0) on the two sides of a forbidden zone. Shifting infinitesimally the point x — iO towards 
the forbidden zone, one can pass it through the cut onto the second sheet. Thus p{x + iO) and p{x — iO) can 
be regarded as two branches of a double-valued function p{x) at the same point on the complex x-plane. One 
can think of e *P^^+*°^ and e'^*^^"**^^ as two independent solutions of the chracteristic (|4.15|) equation for the 
monodromy matrix. Since r2(x + iO) is unimodular, e'^^^"'"''^) e'^*^^"*"' = 1, and the quasi-momentum must 
satsify 

p{x -h «0) + p{x - iO) = 27rnfc, x e Ck (4.18) 

on each of the forbidden zones. This equation represents a Riemann-Hilbert problem equivalent to 1)3. 2|l . if the 
quasi-momentum and the resolvent are related as in (|3.1(l . The integer Uk — rik-i — 1 is the number of (anti)- 
periodic solutions within the fc-th allowed zone, that is the number of times the equation Trfl^x) = 2 is satisfied 
as X moves along the allowed zone between the forbidden zones Ck~i and C^. Equivalently, nk — Uk-i — 1 is the 
number of oscillations cosp{x) makes between the end-points of the forbidden zones Cfc^i and Cfc (see fig. 01. 

We see that, in the classical system, the integrality of B-periods ^A.b\ essentially follows from the unimod- 
ularity of the transfer matrix. Since p{x) is single-valued outside of the cuts, all A-periods are zero. In other 
words, there are no condensate cuts. As we remarked in section|21 it is always possible to choose cuts in such a 
way that condensates "evaporate" and all A-periods of dp identially turn to zero. Such a choice of cuts is unique 
and corresponds precisely to cutting the x-plane along the forbidden zones. The condition that defines the for- 
bidden zones, ImTrO = and |Trr2| > 2, does not necessarily coincide with the condition that — idxDiscp be 
real and positive, which is appropriate for the Bethe ansatz. This is a crucial difference as compared to the 
quantum KdV equation, where Bethe roots condense on the forbidden zones of the classical problem (37) . 

A simple asymptotic analysis of the linear problem H4.13|l . H4.16|l shows that the quasi-momentum behaves 
as p{x) ^ at X ^ 0, this analysis amounts to solving equation 

H4.13|l in the WKB approximation. Hence, G{x) — p{x) + ^ is regular at zero on one of the sheets of the 
Riemann surface H3.3|l . and this singles out the physical sheet of S. 

To complete a proof of equivalence between the Bethe equations and the classical auxiliary problem, we 
should also check that the reality of S(t, a) reproduces the reality condition on the distribution of Bethe roots. 
The reality of S(t, a) implies the following conjugation property of the monodromy matrix: 

W(a;,i,cr) = (T2*(x;i,cr)cr2, Vi{x) ^ (72n{x)a2, (4.19) 

where the bar over matrices denotes complex conjugation of all matrix elements. If Vl{x)i}} — e^^^^^ip^ then 
Vl{x)'4} = e-'P'-''^, and n{x)a2ij = e-'^^'^VaV'- Therefore, aaV- is an eigenvector of fl(x) with an eigenvalue 
Q-iPix)^ i.e. p{x) = p{x), which is indeed equivalent to the reality condition for the resolvent. 

4.2 Finite-gap solutions of the chiral field 

Since the three-sphere is the group manifold of SU{2), the sigma-model H4. 1(1 - 1)4. 3|l can be equivalently 
reformulated as an SU (2) principal chiral field: 

/ Ai + iX2 A3 + iA4 \ / Zi Z2 \ 
9=\ U _ _ ]eSU{2), 

\ -A3 + ZA4 Al~^A2 / V ^1 / 

_ _ _ . (4.20) 

/ ZidZi + Z2dZ2 ZidZ2 - Z2dZi \ 
9 dg = I _ _ _ e su{2). 

\ Z2dZi - ZidZ2 Z2dZ2 + ZidZi 
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The action (|4.1|l takes the form 

(4.21) 



Sam = - / dadr 



where ja — g ^dag = ^Ja'^^ with a = 0,1 ~ t, a, are the right currents. The equations of motion can be 
written in terms of their hght-cone components as 

(4.22) 

d+j- -d-j+ + [j+J-] = 0. 

where the last equaUty follows from the expression for the currents in terms of the matrices g in (|4.20() . 

The sigma-model on 5'^ possesses a global SUlC^) x SU ^(2) symmetry which acts as g{T, a) by left and right 
multiplication by constant group matrices. This symmetry can be identified with a particular 5*0(4) subgroup 
of the 50(6) i?-symmetry group oi M — A SYM. The six scalar fields in the SYM transform under 50(6) in 
the same way as the Cartesian coordinates on the five-sphere in the AdS^ x 5^ geometry. The fields $i and 
$2 have the same quantum numbers under SUlI^) x SUii{2) as the sigma-model target space coordinates Zi 
and Z2 defined in H4.20|) . i.e. ^ Zi and ^2 ^ Z2- Under the left shifts, generated by 

Qt = ^ld<jl^, (4.23) 

where la = gjag^^ = dagg~^, we have g hg, and (^1,-^2) and (^2,-^1) transform as doublets. The 
normalization of the generator is such that Zi and Z2 have Qi = 1- This implies that for an operator 
Tr(*f-^*^ + ...): 

Ql - L. (4.24) 



Under the right shifts, generated by 



Qi-^ldaj^ (4.25) 



47r _ 

g gh, and {Zi, Z2) transforms as a doublet, so that Zi has Qj^ — 1 and Z2 has Qfj — —1. Consequently, 



The scaling dimension of an operator is dual to the energy of the string solution 



3 =L-2J. (4.26) 



27r 



/•2TT 

/ dadrXo^VXn, (4.27) 
Jo 



which is generated by the global time translations. Thus, the Virasoro constraints (|4.3|l become 

lTrjl = lTrf_ = -K\ (4.28) 

These constraints describe a consistent reduction of the 51/(2) chiral field introduced in [7S]. Upon the identi- 
fication 

j+ — inS^ ■ cr, j_ = ik S_ • cr, (4.29) 
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the model reduces to a system of two interacting spins of unit length S+ and S_ 



a+S_ +2kS_ X S+ = 0, 
d-S+ ~ 2kS_ X S+ 0. 



(4.30) 



These equations automatically incorporate the Virasoro constraints. The model has an interesting quantum 
version [791 18UL I5T] . which is a spin system with two spins per site. We do not see an obvious way to relate the 
two-spin chain of [75| to the dilatation operator of A/" = 4 SYM theory (though the Heisenberg model can be 
obtained in a certain limit, see below), but do not exclude the possibility that such a relationship may exist. 

The equations of motion for the principal chiral field can be written as a flatness condition JO . Define a 
current J{x) dependent on a spectral parameter x as 



J±ix) 



3± 



Then the zero-curvature equation 



1 =F a; 

d+J--d-J+ + [J+, J-] = 



(4.31) 
(4.32) 



is equivalent to (|4.22(l . 

On the other hand, it is straightforward to show that H4.32|l is equivalent to the consistency condition of the 
following linear problem 



= dr 



J- 

2 \1~ X 1 + x 



2\l-x l+x 



/ S_i_ • <T S_ • (T 



2 \\-x l+x 
in ( Sj_ • (T S_ • (T 



2 \ 1-x l+x 



* = 0, 

* = 0. 



(4.33) 



The solution to the linear problem now has essential singularities at a; — > ±1, where near the singularity it has a 
solution similar to H4.16|l . and the solutions to the chiral field equations H4.30|l are then defined as in (|4.12|) . We 
are not going to discuss this issue here in detail, since for comparison with section |21 we only need the spectral 
data of the problem (|4.33() . 

Moreover, a more effective way to construct the finite-gap solutions to sigma models was proposed in |84| . 
which is generic for any number of fields in the sigma-model. It uses the functions on the double cover of 
the spectral curve (|3.3fl and its generalizations for larger number of fields. This construction is reviewed in 
Appendix fXl 

The first equation of H4.33|l defines the monodromy matrix: 



ri(x) = Pexp / da 



2tt 



1 



J + 



J- 



2 \ l-x 1 



with the quasi-momentum: 



Tiflix) = 2cosp(a;) 



(4.34) 



(4.35) 



The rest of the story does not differ much from the analysis of the spectral problem for the classical Heisenberg 
model. The only difference is that the Lax pair H4.33|l now has two poles at a; = ±1 instead of one. The standard 
asymptotic analysis yields 

P(^) = -r^ + --- (a^^Tl)- (4.36) 



x±l 
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To express the charges in terms of the spectral data, we expand the quasi-momcntum at zero and at infinity. 
At infinity, C = da + ja/ x + . . and 

Tr ^7 = 2 + ^ da,da2 Tr ^0(^1)^0(^2) + . . . = 2 - + ... = 2- ^^'^^ ~ ^"^^^ + . . . , (4.37) 

Here we assume that the classical solutions describe highest-weight SUl{2) x SUii{2) states. For highest weights 
in large representations, \Ql,r\'^ can be replaced with (Q|)^ = and {Q%)^ = [L — 2J)^. Thus we have 

, , 27r(i- 2J) , , , ^ 

p{x) = + ... (x ^ 00). (4.38) 

V A X 

At X ^ 0, £ = So- + ji — xjQ + . . ., which can be written as £ = g^^{da — xIq + . . .)g. Then, 

n{x) = .g-^(27r)Pexp (^--^^ dalo + ..^ .g(0). 

Because of the periodicity of g{cr), ri(0) — 1 and p(0) — 2nm. Expanding further, we get 

Tvn = 2 + — daida2Tvlo{ai)loia2) + ... = 2 ^ + . . . = 2 + . . . . (4.39) 

2 Jo A A 

Hence, 

27rL 

p{x) =2T:m+ —^x + (x ^ 0). (4.40) 
V A 

By the same arguments as in the previous section, the quasi-momentum is a meromorphic function on the 
complex plane with cuts, satisfying equation 13.2|l on each cut. Subtracting the singularities at a; — > ±1, we get 
the resolvent 

Gix)=pix) + 7 + — T, 4.41 

a; — 1 a; + 1 

Because the poles at a; = ±1 cancel, the resolvent is an analytic function on the physical sheet and, as in H2.12|l . 
it can be represented as an integral of a positive density 



G{x) =Jd^f^^. (4.42) 



The proof of this spectral representation is the strandard argument based on analyticity of G{x). We define 
2iTTp{x) = G{x + iO) — G(x — iO). Then the right-hand side can be represented by a contour integral with the 
contour surrounding all the cuts. The only singularity of the integrand on the outside of the contour is a pole 
at y = X with residue G{x). Shrinking the contour, we get 14.42|l . 

The asymptotic behavior of the resolvent, G{x) ^ 2t:[k — (L — 2 J)/%/A]/a;, translates into the normalization 
condition for the density: 

27r 

dxp{x) = ^{A + 2J - L). (4.43) 



dx^^ 27rm (4.44) 



The density is subject to two extra constraints, which follow from the asymptotic behavior of the resolvent at 
zero: 

_±IG{x)dx _ f^^pix) 
27ri 
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and 



1 



G{x)dx 



— dx 



p{x) 



2tt 

7x 



(A-L). 



The quasi-momentum satisfies (|4.18|) : 

p{x + iO) + p{x — iO) = 27rrifc, 



xeCk 



(4.45) 



(4.46) 



This equation is a simple consequence of unimodularity of the transfer matrix. Taking into account the spectral 
representation H4.42|l . it can be recast into an integral equation for the density: 



G{x + iO) + G{x - iO) = 



2ttk 



2ttk 



27mj, 



X e Cfe. 



(4.47) 



We obtained again the same Rieniann-Hilbert problem as for the long spin chain (|3.2|l or for the classical 
Hisenberg ferromagnetic (|4.18|l . but with a different pole-structure. This equation is a direct generalization of 
H2.15|l and can be called the classical Bethe equation for the chiral field. It can be solved in the same way as 
H2.15|l . Namely, we look for a differential dp defined on the hyperelliptic surface having double poles: 



dp ~ dx 



the eq. 



at X — ±1, and behaving as dp 



dp ~ TTK 



dx 



y+ 



^(L-2J) 



y- 



dx 



at X 



oo. We write, generalising 



(.T-i)2 {x + iy 



1 X+l 



bkx 



fe-1 



(4.48) 



fe=i 



where y± 



y\x=±i^ y± — dx 



x=±l 



and the coefficients bk are determined, as in H3.6|l by vanishing of A-periods. 



The challenge is to guess what quantum Bethe equations can reproduce (|4.47|l and its solution (|4.48(l in the 
scaling limit. This does not look like an easy problem. The Bethe equations for the Faddeev-Reshetikhin model 
|79| is an obvious candidate, and not surprisingly - this model was constructed as a quantization of the chiral 
field. The Bethe equations of 79- indeed reduce to 1)4.47(1 in the scaling limit. The parameters, however, do not 
match literally. 



4.3 Comparison of string theory to perturbative gauge theory 

We are now in a position to compare the results of the sect. [21 for the spin chain, that describes the one- and 
two- loop perturbation theory of the JV — A planar SYM gauge theory, with the sigma model of the dual string 
theory. We should compare either their generals solutions (|3.27|l and (|4.48l) , or the general equations (|3.21() and 
(|4.47|l . together with the asymptotic conditions for the resolvent G{x) or quasi-momentum p(a;). 

First we will try to get a qualitative idea at one loop, by comparing the small X/L limit of sigma model to 
the classical magnet introduced in this section as a naive continuous limit of the XXX spin chain. 

The general solution of the chiral field should reproduce the solution of the spin chain in the weak-coupling 
limit. On the classical level one can take the limit directly. If the 't Hooft coupling is small, then n = A/VA is 
large and we can expand the solution of (|4.30() in the inverse powers of this parameter: 

S.-S±isxS.-Ji,S + 0(-i). (4.49) 
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The 0(1/k^) term ensures the correct normalization of the classical spins: Sj_ = 1. The world-sheet time 
variable should also be rescaled such that dr — *■ 9r/(4K). Then 

d± = ■^dr±da, 

and the leading 0(1) term in the equations (|4.30|l cancels identically. The next 0{1/k) term gives the equation 
of the Heisenberg magnet: 

Sr^Sx S„„. (4.50) 

This procedure is in fact equivalent to the derivation of the Heisenberg equations from the sigma-model of 
ref. which is based on separation of fast and slow modes of the string and also involves rescaling of the 
world-sheet time 

Plugging (|4.49l) into H4.33I) . rescaling x — * Attkx and taking k to infinity 



iS cr fl 

c^d.- - — + o - 

AlTX \ K 



/ ^ (4-51) 

iS ■ cr i[S X So-) • cr ^ / 1 



AkM ^dr- - + O , 

'm'^x'^ Attx \k 

we recover the Lax pair of the magnetic. It means that on level of classical integrability, some finite-gap solutions 
of sigma-model on string side turn into the classical solutions of the Heisenberg spin chain. 

However, to make quantitative comparison between the prediction of SYM and string theory for the anoma- 
lous dimensions (|3.25|l and H4.45II up to two loops, we have to compare the solutions to basic integral equations 
(|3.21|l and (|4.47|l . together with the normalizations of densities and zero-momentum conditions. In other words, 
we have to see how the geometric data depends on the integrals of motion. In order to do this, let us take 
directly the weak-coupling limit in the integral equation H4.47|l . Rescaling for example a; by a factor of 47ri/\/A, 
the normalization conditions 14.43|l . (|4.44|l and l|4.45(l become 

, , , A J1-J2 J , A-L 
ax p(x) 



2L 2L L 2L ' 

G{x) = fdx^^ 27rm, (4.52) 



2TTix J X 



and the integral equation (|4.47|) transforms into 



+27rnfc, x&Ck. (4.53) 



^^We are grateful to A.Tseytlin for discussions of these issues and the subtleties concerning the general two-loop comparison, as 
well as sharing with us unpublished results on a similar approach to relating sigma-models with classical spin systems. 
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In the limit A/i^ ^ 0, A = i + 0{X/L), and these formulas exactly reproduce the scaling limit of Bethe equa- 
tions (|2.15|l for the one loop perturbative SYM, together with the relations (|2.13I) . H2.16|I and H2.18|l completely 
defining the anomalous dimension. We see that the one-loop physics (at least in this SU{2) sector) of the SYM 
is immediately reproduced from the string sigma-model for very general solutions. 

Comparison to two loops takes a little more work. If we expand the rhs of H4.53II in powers of T = le-i^iL'i ' 
see that the coefficient of the term is half that of the corresponding coefficient in H3.21|l . This suggests that 
one should replace the spectral parameter in (|3.21|) with x — x + T/x + 0{T'^). The string density is thus related 
to the gauge theory density by a change of variables, Pg{x) — p{x — T/x + ...), where the subscript g refers to 
the gauge theory solution in section 3.2. To reproduce the results of section 3.2, we need to change variables in 
the integral equation and in the normalization conditions. The change of variables in the momentum condition 
in (|4.52|l gives 

Jdx^(^l + ^ + ..}j^2nm, (4.54) 
which is the same as (|3.23|l . The energy conditions transforms to 

A-L = 2TL J dxpg{x)(^-^ + ^ + , (4.55) 

in agreement with ()3.25f) . The normalization of the density becomes 

, . f T\ J A-L , , 

dxpg{x)[l + ^j=- + ^^. (4.56) 

The unwanted terms cancel in virtue of H4.55|l in agreement with the canonical normalization of the gauge 
theory density (|2.14l) . Finally, 



f dypjy) ^ f dz pgjz) f Pg{z) _^ 

J X - ^ - y J X - z x J z"^ 



Using the energy condition (|4.55() . we find that Pg{x) satisfies H3.21|l with the two-loop accuracy. We see that 
the sigma model matches the gauge theory prediction up to two loops for all long XXX operators. The gauge 
theory resolvent is related to the resovent in the sigma model in a simple way. At one loop they are the same. 
At two loops, an 0{T) correction arises because of the change of variables: 

G,(x) = G(x)-r^^M-^-t-... . (4.57) 
This equation describes the map of higher charges in the sigma model to those in the gauge theory. 



5 Examples 

Let us now turn to particular examples of our generic solution. We consider first the BMN case, which from 
integrable system point of view corresponds to a solitonic degeneration of our general finite-gap solution. Then 
we discuss in detail the solutions with one and two cuts. 
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5.1 BMN states 




Let us show how the BMN formula for the anomalous dimension jlOj is reproduced from the integral equation 
H4.47|l . The BMN limit corresponds to a vanishingly small macroscopic density p{x), except for contributions 
concentrated around the zeros of the r.h.s. of H4.47() : 

(5.1) 

In terms of the filling fractions (|3.1U|I . the normalizations (|4.43|) . 14.4411 and (|4.45|) in this approximation become 

i:s..^|(A+2.-.). x:|.o, i:|^^(A-i)^ ,5., 

Let us introduce the quantities Nk, which parameterize the filling fractions as 

(5.3) 

The Nk 's have the meaning of occupation numbers of string oscillators and should be large for the semiclassical 
approximation to work. On the other hand the filling fractions Sk are small in the approximation we use, so the 
coupling must satisfy a/A ^ Nk- The second of the equations (|5.2(l then leads to the zero- momentum condition 
in the form 

J2Nknk = 0. (5.4) 

k 

The total momentum vanishes in this case, since all Sk are negligibly small, which is consistent only for to = 0. 
The sum and difference of the first and third equations of l|5.2|l gives rise to 

(5.5) 





which are the familiar BMN results, except that L is replaced by A, but which only leads to corrections beyond 
the BMN limit. From the point of view of integrable systems the BMN limit corresponds to the solitonic 
degeneration of the finite-gap solutions. In other words, the cuts shrink to points. One can also consider partial 
degenerations when some filling fractions are small but others are arbitrary, which corresponds to solitons on 
the background of finite-gap solutions. These solutions describe small oscillations around macroscopic rotating 
strings. 



5.2 Rational Solutions 

In this subsection we will consider solutions with a single cut. These these solutions are dual to a class of 
semiclassical string solutions described in l^^. We will demonstrate that this simple example is not only 
instructive for our general approach, but also show that all rational solutions have two-loop agreement between 
the SYM prediction and the semiclassical string prediction. 
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= 2TTn + 2-f ^^^^ =27rn + G(a; + iO) + G(x-iO), x€C (5.6) 
J c 2; — ^ 



Spin chain solution for the one loop gauge theory 

First, we consider the solution to the quasiclassical Bethe equations of section 13 With only the one cut C, 
eq. I|2.15|l reduces to 

c 

The momentum constraint (j3.12|l imposes the further condition that 

na — m (5.7) 

where m is an integer. With a single cut, it is straightforward to find the resolvent G{x) from H5.6|l . where one 
finds 

G(x) = i f - - 27m + -J(2TTnx - + Smnx] (5.8) 
2 \x X J 

i.e. the curve in this case is rational, and can be defined (for n 7^ 0) by the equation 

y' = ^ ((2-- - 1)^ + 8^) - - i (1 - v) ^ + ^'-'^ 

The coefficients inside the square root are chosen to eliminate the pole at a; = and to reproduce the large x 
behavior of G (EISJ), 

G{x) ~ ^ + ... (5.10) 

X — >oc X 

With the resolvent in hand, we can now find the anomalous dimension, which is proportional to the linear 
coefhcient of the Taylor expansion (|3.13(l of G{x) about a; = 0. Hence, 

A / dG\ \m(n~m) 

Comparing this solution to the string solutions in [^Hl, we see that this is the same value as the circular string 
solutions with J3 = and Ji and J2 satisfying 

Ji + J2 = i 

(5.12) 

mJi + (m — n) J2 = 

Notice further that the special case of n = 2 and m = 1 is the original solution of [20]. In other words, the 
Frolov-Tseytlin solution may be thought of as a BMN type solution where there are creation operators a^^j 
acting on the chiral primary ground state. From this, it is intuitively clear why this state has a higher energy 
than the folded string with the same i?-charges; the folded string is made up of -j creation operators a^^^ and 
■J creation operators aLi- 

The number of Bethe roots is limited to half the number of sites in the chain. Therefore, a — m/n < ^. 
However, the resolvent (I5.8|l is well behaved if m > n/2 and 7 is invariant under m ^ n — m. This suggests 
that there is another configuration of Bethe roots which leads to the same resolvent, but with a for the new 
configuration given by 1 — a of the old configuration. It is easy to see what the new configuration is by considering 
the position of the branch points of H5.8(l which are located at 



1 / n — 2m 



x+ — — 
^ 2tt 



± — \/ — nm | (5.13) 
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If m < n/2 then the real part in (|5.13|l is positive, but it is negative if to > 7i/2. In fact, it is clear that under 
the transformation n —^ —n, m —n + m, G{x) is transformed as G{x) — > —G{—x) (this transformation leaves 
the curve (|5.9|) intact). Hence, if we start with a configuration with m > n/2, this is equivalent to having 
n' = —n, to' = —n + m and a — m! jn' < i. In other words, this corresponds to BMN states with L{n — 'm)/n 
creation operators aL„ acting on the chiral primary state. 

In the case of n = 2to, the solution considered here is similar to the 50(6) singlet solution constructed in 
|24| . But we could have also chosen the Bethe roots to consist of a condensate between the branch points and 
two tails leading to ±ioo, this is the solution originally constructed in [22]. To see how this works, recall that 
p{x)dx has to be positive definite along the distribution of the roots. The special case n = 2m has branch points 
at x± = ±i/(27m) and locally from the branch points p{x) ~ ^/x± — x. Hence there are three directions coming 
out of the branch point on which the roots can lie. Two of these directions correspond to a cut connecting 
the two branch points on either side of the imaginary axis. The third direction is along the imaginary axis out 
to infinity. This would normally lead to a nonnormalizable density. However, we can also choose the opposite 
sign of the square root and add a constant density along the imaginary axis. This then leads to a condensate 
between the branch points and two normalizable tails extending out to infinity. 

Solution of the classical magnet 

The simplest rational solution of the classical ferromagnet equations in (|4.4() is 



WT + na, ~ 00, with 



w = cos ( 



(5.14) 



For this solution the total angular momentum is 



L f^^ L 
Sz = — j da COS0 = — cos 00 (5.15) 
47r Jq 2 



and the anomalous dimension or Hamiltonian is given by 



A „ A r^"" , ..9 . 9 ,.9n An2 



7 = - ^ da [id,0f + sin2 0{d,^f] = ^ sin^ 9o (5.16) 

Now, for this particular solution one has for U defined in H4.1(J|I 

/ cos6io sin 6^0 e""" \ 

^l*=o - . , , (5.17) 

\ sm^oe -cos 6*0 / 

and the first of the equations of the auxiliary linear problem in (|4.9|l can be easily integrated: 

«'(x;0,a) = -— 1— 2:. V J ^ ^^^^^^ 

\ 2x \ J ^ 
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where 



k± = Tm 



cos 9o 1 
- 1 ± - 

2'Knx X 



cos e/o 



1 



p± — im \ ±— \ x^ — 

\ X \ Tm 



(5.19) 



cos Wo 



1 



Att^t 



are functions on H5.9|l . Putting cr = 27r, taking the trace and observing that cosp^ — cosp_, we get 

Tr*(x;0,27r) = Trf7 = 2cosp± (5.20) 

Hence, the quasi- momentum on the physical sheet is p+ in H5.19|l . 

We now have to impose the momentum constraint. One way to find the value of the momentum is to 
compute the expectation value of a shift operator for a quantum state. The correspondance principle then 
relates this to the classical value. The quantum state is a collection of L "up" spins, where "up" for the j*'' 
spin is with respect to a polar and an azimuthal angle 9j and . Hence the state has the form 



I state) 



v 



y / 



cos — cos — 



2 / 



V 



(5.21) 



and it is parameterized, as usual, by spinor "half-angles" 6/2 which in our case correspond to the sigma- model 
co-ordinates given below. 

The shift operator e*^ sends the spin at site j to the spin at site j -t- 1. Hence, the expectation value of the 
shift operator is 



(state|e^^|state) ^ [] ( cos | cos ^ + sin | sin Z^eK*^+i-*^) 



(5.22) 



]J h + i sin^ ^ ) » exp i ^ sin^ -j-A(b, 



where we have assumed that Adj << 1, A(/)j << 1, which is valid in the long wavelength limit. Hence, the 
classical momentum is given by 



P 



dcrsin^ —drr<t> — 27rnsin'^ — . 
2 2 



The momentum is required to be a multiple of 27r, which puts the restriction on 

2m 

cos t'o = 1 

n 

showing that p+ in (|5.19|l is consistent with G{x) in l|5.8|l . 

Let us now write everything in terms of the _R-charges Ji and J2 . We have that 

Ji = ^ + 53, J2 = \-S^ = ^{1 - cos^o) = Lsm^ ^. 



(5.23) 



(5.24) 



(5.25) 
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We also have that a = hence, P is given by 

P ^ (5.26) 
and the momentum condition then forces an = m where m is an integer. This then leads to 

which is the same as ()5.11(l . Note that the singly wound Frolov-Tseytlin circular string reduces to a doubly 
wound classical solution of the XXX ferromagnet. 

Solution of the sigma-model 

The string is restricted to live on S^, so we choose the parameterization 

Zi=Xi+ 1X2, Z2 = X3 + 1X4, (5.28) 

where 

|^i|' + |^2|' = l. (5.29) 

Let us consider solutions of the form 

which, in particular, corresponds to the rational case (jA.lOjl of the general construction of Appendix The 
angle Oq/2 is the polar angle on S3 in generalized coordinates and so Oq ranges between and tt. The equations 
of motion for Z12 lead to the relation for the parameters of the solution 

2222 /'c;01^ 

Wi — nil = ■W2 — (5.dl) 

while the Virasoro constraints H4.3|l require 



= {wl + ml) cos^ — + (i«2 + ml) sin^ y 

2 ^0 , . 2 ^0 n 

wimi cos — + W2m2 sm — 0. 



(5.32) 



in particular meaning that 



and 



^^^2 ^ ^ W2m2 ^.^^2 ^ ^ ^irn, 

2 W2m2 — wimi 2 wimi — W2TO2 

^2 _ (W1W2 + miTO2)(wim2 - W2TO1) 
W2'm2 — wimi 



which is, together with (|5.31|) solved by 

1 



wi — — \ \/Wi 



2 V VWi^ 

m/ = - , \/Wi 

2 \VWJ ^ 



(5.35) 
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which is equivalent to (|A.11|) . Using (|5.35|) . the formulae in (|5.33ll become equivalent to ljA.14|l . 

The energy, which is A, and the i?-charges of the string are given by 



_cos--u;i=yAcos'y«;i ^^gg^ 



Ji = \/A / 
Jo 

J2 = \/A / ^ sin^ = %/A 

Jo ^'^ ^ 

Now, the relations H5.31|l and H5.32|) lead to the equations 

{w\ + ml~ ml) ^Ji — "v/Awi^ — {J2Wi)'^ — 



• 2 ^0 

Sm yW2 



(5.37) 



JlTOi + J2m2 = 

Note that the situation simplifies greatly if mi — — m2 = m, in which case we find that 

A^ = (Ji + J2)2 + Ato^. (5.38) 

In order to compare this to a single cut solution of the Bethe equations, it is convenient to rescale x by a factor 
of a = Equations (|4.52(l are then modified to 

'^^^(^) = A + ^^' 

G{x) ^ 27rm, (5.39) 



27rja; 

-8;ffe/2^^(^) = ^-^' 

and (|4.47|) becomes 

G(x + iO) + G{x - iO) = 2-1 d£, + - 27rn, (5.40) 

J X — ^ 2 \x — a X + a J 

Then the general form of G{x) for mi = m, m2 — m — n, satisfying H5.40|) and compatible with the general 
solution 14.4811 for the chiral field, is 



Gix) 4 + ^) + i f^i^^ + ^^^^ ) V Ax^ + Bx + G -nn, (5.41) 

4: \x — a X + a J 4\ x — a x + a y 

where, in order to cancel the poles at a; = ±a on the physical sheet, we must satisfy the relations 

B = Anne 
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In order to have the asymptotic behavior of 12.13|) we also require 



and to satisfy the momentum condition (|4.44|) we have 



VA ( + ^ ] = 47m, (5.43) 



G(0) = ttkVC ( ^^== - ] -Tm = -27rTO. (5.44) 



If we define Y = \/l — e^, then equations (|5.42|) , H5.43|l and (|5.44|) lead to the equation 



2 \1 + Y 1-Y 

Using and the definit ion for Y, one can show that asymptotically, G{x) in H5.41|l behaves as 



a;— ^oo 



Substituting H5.45(l for and comparing H5.46|l with H5.39fl leads to the equation 

n2(l - Yf - (n - 2m)2(l + Yf = ^I^i^ll—^^T^. (5.47) 

vA 

It is then a straightforward excercise to show that H5.39(l is consistent with H5.45|l and l|5.47(l if 

Ji = J, = "^L. (5.48) 

n n 

Hence, we have that m ~ mi, n ~ nii + m2 and Y satisfies the quartic equation 

(^2(1 _ Yf - {n - 2mf{l + Yff - 16(n - 2m)2(l - Y"^)^ = 0. (5.49) 

A 

In the special case where n = 2m, then Y — 1 and H5.45|l reduces to H5.38|l . We also see using (|5.38l) that G{x) 
in (|5.41|l simplifies to 



Finally, expanding G{x) about x ~ and matching the linear term to the negative of the third equation in 
H5.39|l . one immediately finds |E 



Solution of the spin chain at two loops and comparison with the sigma-model 

For general values of n and to, we can at least solve for A in a series expansion of X/L^ . Solving for Y in 
we find 

(n-2TO)2 A 3(n-2TO)4A2 {6{n ~ 2mf - n^){n - 2771)"^ 
2 I2+ 8 l4 16 
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Substituting this into 15.45|) wc find 

m(n — m) A m{n — m)(n^ — 3m(n — m)) 
2 I 8 l3 



(5.52) 



The term hnear in 7 in (|5.52|l is clearly consistent with l|5.11() but now let us consider the two loop term. To 
do this, we consider the one cut solutions to the two loop approximation of the Inozemtsev chain. Using H3.21|l . 
we can write down a general solution for G{x) compatible with (|3.27|l as 

1 T f D E T 



G{x) = ^ + ^+ ( — + — + — ) ^^2;2 + Bx + l- Trn, (5.53) 
where we recall that T — iQ^i^^i ■ Cancellation of the double and single pole at a; = leads to the equations 

£;+iBr = 0, (5.54) 

D + EB + AT --B^T =1. (5.55) 
4 

The momentum condition (|3.23|) to linear order in A/L^ gives 

\dB + -EB'^ - -EA - —TB^ + -TAB - T ( —B^ - -AB] ^ 7r(n - 2m), (5.56) 



4 8 2 16 4 V16 4 , 

where the parentheses term comes from the second term in H3.23|l . Finally, the constraint on the constant term 
at infinity is 

D\^=2Trn. (5.57) 
Solving these equations up to terms linear in \/L^, wc find 

A 



D = l-S, E=^-^{n-2m), A ^ {27mf {1 + 2S) , 
B = -47r(ri - 2m) (l + S + 4m(?i - m)-^ 



(5.58) 



where 

(5 = -— -j(n^ - 6m(n - m)). (5.59) 
Inserting these values into (|5.53|) . we find that 

G'(0) — — 47r^m(n — m) + 27r^m(n — m){n — 3m)(2n — 3m) — ^ + ... 



— G"'(0) — — 167r^m(7i — m){ii? — 5m{n — m)) + 



(5.60) 



Applying this to (|3.25() . we finally obtain 

A ^//„N 3A2 G"'(0) m,{n — m) X m,(n — m)(n2 — 3m,(n — m)) A^ ('^p,^ ] 

^"~8^^^" 128^^3 3! 2 L 8 L3 + • • ■ ' > 

which, finally, agrees with (|5.52|) . 
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5.3 Pulsating Solutions 



This section is somewhat outside the main development of this paper, but it indicates that similar techniques 
may be applied to operators that are not of the XXX type. The relevant observation is that the isometry 
group of the sigma model on is SU{2)ii x SU{2)l, while the isometry group of the Heisenberg magnet is 
SU{2)ii X U{1)l- Hence, there will be classical solutions of the sigma model that do not reduce to solutions of 
the Heisenberg magnet as the coupling is taken to zero. The solutions that do reduce to the magnet are those 
where Ql equals the bare dimension of the operator. 

We now relax this condition and consider the so called pulsating string solutions. These string solutions are 
not dual to XXX operators, but instead are dual to operators comprised of all six types of scalar fields. To 
proceed we consider 

Zi = cosxe*""", Z2 = sinxe^"^, (5.62) 

where x aud (j) are functions only of r. The Virasoro constraints 14.3|l give 

= X^ + rn^ cos^ X + 0^ sin^ x- (5.63) 
We also have from the equations of motion that 

—X tan X ~ + = X cot X — — <j>'^ + '^^x4> cot x + i0. (5.64) 
Setting the imaginary piece to zero and integrating leads to the equation 

■ j/Vx 



■ 2 ' 

sm X 



(5.65) 



where the integration constant is determined by the _R-charge J. If we now insert (|5.62|l into 14.20|l and (|4.33|) 
we find the auxiliary linear equation 



d„^ix-0,a) = ^^{ : ^ )*(x;0,a), (5.66) 

~ 1 \ -e"""x iJ/^ 



where we have assumed that 4'\r=o ~ cosxI^-^q ~ Integrating this equation and using the Virasoro constraint 
in (|5.63|l we find that the quasi-momentum is given by 

p{x) — — -]j ^"^(^^ ~ 1) ^ VA^) + — ^ — TTTO . (5.67) 

and the resolvent G{x) is 



G{x) = [kx+J f ^(a;2 - 1) - -^x] + ( - ^ ] x^ I - mn . (5.68) 



1 I V V2 ' ' J V ^ 

Rescaling again x by Attk, we get for G{x 




2 



= a: + wf2^m (x^-—^]-{l^f3)x] + P{2 - P)x^ \ - nm . (5.69) 
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where /3 = 1 — ^. When A ^ this becomes 

Go{x) = ^ + V (27^"!^ - (1 - Pyf + P{2 ~ P)^ - irm . (5.70) 

In [25, it was shown that the resolvent coming from the one- loop Bethe equations for [0, J, 0] representations of 
S0{6) was proportional to Go{x) — Go{—x). In terms of (|5.7U|) we see that this corresponds to adding together 
the resolvents for two classical solutions with ±m. This suggests that the full resolvent, dual to these states, 
is G{x) — G{—x) where G{x) is given by (|5.69|) . It is also interesting to point out, that at the limit A ^ the 
elliptic curve of H5.69|l turns into the rational of (|5.7U|) , which is reminiscent of a corresponding phenomenon in 
softly broken A/" = 4 Seiberg-Witten theory. 



5.4 Elliptic Solutions 



The next case in simplicity after the rational solutions of Bethe equations are solutions with roots distributed 
on two symmetric cuts. The Riemann surface (|3.3|) of a two-cut solution 



y 



I 2 2\/ 2 2\ 

[x -Xi)(a; -X2) 



(5.71) 



is a torus and is just a particular K = 2 case of (j3.16|l . The resolvent for the two-cut solutions can be expressed 
in terms of elliptic integrals. Elliptic solutions are parameterized by a pair of integers, n and m, which determine 
two independent periods of the differential dp. Two types of elliptic solutions found in |32| are dual to circular 
and to folded strings rotating in . We will show these solutions essentially exhaust all possiblities. In the 
process we will show that not all pairs of integers n and m are compatible with the reality condition and that 
some pairs are related by symmetries and correspond to one and the same distribution of Bethe roots. 

According to sect.|21 two-cut solutions of Bethe equations are constructed with the help of a mcromorphic 
differential on a double cover of the complex plane with two cuts (xi, X2) and (— X2, — xi). The differential with 
the desired properties is 



dp — 



dx{\ 



a - 



X1X2 ' 



(5.72) 



v/(x^-a:2)(x2-x^) 

The constant term in the numerator assures the correct normalization at infinity (|3.11l) . The square root 
is defined in the standard way, such that it coincides with its algebraic value at infinity. The value of 
(xp- — x^)(a;^ — x|) at a; = depends on how the complex plane is cut. We assume that ■\/ (0^ — x^ ) (0^ — x^ ) = 
— X1X2, which is the standard convention for real positive xi and X2. The anomalous dimension H3.13|l in this 
case is 

^ - 2a X? + x^ ^ 



7 



Wtt^L \ X1X2 



2x^X2 



(5.73) 



The periods of dp can be expressed through the complete elliptic integrals of the first and second kind. In order 
to convert them into the standard Legendre form, we first compute the integrals for positive real xi and X2 and 
then analytically continue to complex values of the parameters: 



dp = —2i 



dp^ 4 



dx ■ 



X1X2 

2x2 



dx ■ 



- x'^)(x^ 
a 



X?) 



X2 



{l-2a)K{l-r) - 



E{1 - r) 



X1X2 



\J [xP- — xf )(a;^ — X2) 



2 

X2 



(1 - 2a)K(r) + 



E{r) - K{r) 



(5.74) 
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where the modulus of cUiptic integrals is 

r = ^ (5.75) 

Our conventions about elliptic integrals and some of their useful properties are listed in the Appendix ^ 

If we assign the momentum quantum numbers n and ~n to the two cuts, so that the B-period of dG is 
normalized to Attu, and put a condensate between the cuts, such that the A-period is equal to 2Tnn. The 
conditions (|3.14fl then become 

Ed - r) 

(1 - 2a)K(l - r) ^ ^ ' = 27rimx2 

^ (5.76) 

E(r)-K(r) 

Vr 

Exluding X2 from these equations we find 

^_2a=l_ nE{l-r) + im{E{r)-K{r)) 
^/r nK{l — r) — imK{r) 

This equation determines the modular parameter r as a function of the filling fraction a. Once r is known, we 
can express the end-points of the cuts, and therefore the resolvent and the anomalous dimension, as implicit 
functions of a. Using the Legendre relation (jB.3|) in intermediate calculations, we find 



1 1 

Xl 



4 nK(l - r) - imK(r) 

^ \ (5.78) 

^ nK{l — r) — imK{r) 

The above equations in principle solve the problem of computing the anomalous dimension as a function of the 
filling fraction. However, the anomalous dimension must be real which is not always true for arbitrary n and 
m. The reality condition constraints allowed values of n and m. There are also certain symmetries that relate 
solutions with different n and m. 

The modular parameter r is an auxiliary variable in the above equations, it just parameterizes the dependence 
of xi and X2 on a. Nothing changes if r is replaced with 1/r. The elliptic integrals tansform linearly under 
such transformation, according to ljB.21) . As a consequence, (|5.77|l . (|5.78|) have almost an identical form when 
expressed in terms of 1/r. A short calculation shows that the replacement of r by 1/r results in a simple 
transformations of n, m, xi and X2: 

n — > n — m, 

Xl ^ X2, (5.79) 

X2 Xl. 

The resolvent and the anomalous dimension are manifestly symmetric under interchange of xi and X2. Therefore 
two pairs of integers, (n, m) and {n — m, m) correspond to one and the same solution of the Bethe equations. 
Hence, n is only defined modulo m. 

Complete elliptic integrals K{r) and E{r) are multivalued functions of r. They have a logarithmic branch 
point at r = 1. We can pick a different branch of K{\ — r) and E{1 — r) by pulling r around zero. K{r) and E{r) 
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are analytic at zero and do not change, but K{1 — r) and E{1 — r) pick up discontinuities at the logarithmic 
cut. The discontinuities can be easily calculated from the defining integrals (|B.1|) : 



K{1 -r-iO)- K{1 -r + iO) = -2iK{r) 

and 

£;(! - r - iQ) -E{l-r + iO) = 2i {E{r) - K{r)) . 

As a result, 

K{1 - r) ^ K{1 - r) - 2iK{r), 

(5.80) 

E{1 - r) ^ E{1 ~r) + 2i {E{r) - K{r)) . 

This transformation, when applyed to (|5.77|) . H5.78|l . shifts m by 2n. This is not the only effect of the analytic 
continuation. The ^/r also changes sign when r encircles zero. Therefore X2 and 1 — 2a will change sign: 

m — > m + 2?!, 

Xl Xl, 

(5.81) 

a I — a. 

To get the physical filling fraction (a < i) after the transformation, the initial solution should have a > ^. It 
thus corresponds to an analytic continuation of some eligible solution beyond a = ^. The anomalous dimension 
(|5.73|l is manifestly invariant under (|5.81l) . The resolvent is also invariant because the change of sign of the 
numerator in l|5.72|) is compensated by switching to another branch of the square root in the denominator. 
Indeed, the cuts connect xi with — X2 and X2 with — xi after the transformation, rather than xi with X2 and 
— X2 with — Xl. Such a rearrangement of the cuts changes the sign of the square root {x^ ~ x^){x'^ — x|). The 
symmetries (I5.79|l and ()5.81(l correspond to transformations in the homology basis: {A, B) (A, B — 2A) and 
{A^B) {A + B,B). While the former transformation is a genuine symmetry, the latter changes the filling 
fraction and maps a physical solution of Bethe equations to an analytic continuation of that or of some other 
solution. 

Apart from integrality of periods of dp, solutions of Bethe equations must satisfy the reality condition on 
Xl and X2. There are four possibilities: (i) xi — 5c2, (ii) xi, X2 are pure imaginary, (iii) xi — — X2, and (iv) xi, 
X2 are real. There are no solutions with real cuts, which excludes case (iv), but it is known that such solutions 
exist for spin s — Heisenberg chain that is also relevant for the AdS/CFT correspondence and describes 
strings rotating in and AdS^ simultaneously |16| . The other three cases are possible for various n and m. 
Careful inspection of eqs. H5.77|l . H5.78|l shows that there are three series of solutions that satisfy the reality 
condition (see Appendix lUl for details): (i) m = 0, (ii) n = and (iii) m = 2n. The first two cases were studied 
in and it is straightforward to check that the general formulae here reduce to those in when either m or 
n is set to zeio^^. We should mention that only the minimal case of to = 2 has been discussed in [3111111 ■ The 
reason why Bethe state with to = 1, which has lower energy, is not minimal is that it violates the momentum 
constraint, since the total momentum of Bethe states in case (ii) is ttto. 

The case (iii) can be obtained by applying tranformation H5.81|) to case (i) which should be analytically 
continued past a = ^, because the transformation maps a to 1 — a. It was conjectured in |33| that such an 
analytic continuation describes an "umbrella" -like distribution of Bethe roots, with the condensate along the 

^^It is easier to compare to j46|. where all elliptic integrals are reduced to the standard Legendre form. 
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imaginary axis and two arrays of roots at both ends of the condensate which are approximately parallel to the 
real axis (that is, with xi — — X2). Such states were found numerically for small lattices The solution with 
m = 2n has all the expected properties: it is equivalent to an analytic continuation of the solution with m = 0, 
it contains a condensate, and the tails of the condensate have end-points xi = — X2. We thus conclude that 
the elliptic solution of Bethe equations with m = 2n is the "umbrella" state conjectured in It was also 

observed in \32\ \'6'6\ that the solution with m = (case (ii)) is symmetric under a — s- 1 — a. From the point of 
view of the general approach developed here, this is a consequence of the invariance under the transformation 

6 Discussion 

In this paper we have discussed a general, integrability based approach for connecting the two sides of the 
gauge-string duality. On the gauge side we have composite operators, whose renormalization is governed by a 
quantum integrable spin chain, while on the string side we have a classically integrable sigma-model. 

Our explicit calculations were mostly restricted to the SU{2) subsectors of SYM and string theory. The one- 
and two-loop operator mixing in the SU(2) sector can be resolved for arbitrarily large operators and expressed 
in terms of simple geometric data. A particular operator is characterized by a hyperelliptic complex curve, 
endowed with extra meromorphic differential, whose periods obey certain integrality conditions. It turns out 
that solutions of the sigma-model are expressed in terms of the same analytic data, with a slight difference 
in the pole structure of the differential. This difference however accounts for all orders in the BMN coupling 
X/L^. On the gauge theory side, the geometric data is derived by solving the scaling limit of the discrete Bethe 
equations, which are themselves exact, and not restricted to operators of large size. The discrete counterpart of 
the geometric data in the sigma model, if it exists, would describe the full non-perturbative spectrum of single- 
trace operators in the SYM theory. Presumably, it can be derived from quantum string theory in AdS^ x S^, 
about which little is known. Our results suggest that the quantum sigma-model can likely be solved using a 
Bethe ansatz whose scaling limit produces the integral equation (|4.47|) . As sometimes happens for integrable 
systems, the full quantum solution of the model may even be uniquely determined by its semiclassical limit and 
an underlying analytic structure. 

The energy of the string solutions is analytic in the 't Hooft coupling A and can be expanded in a power 
series in X/L^. This expansion, though being very nontrivial on the string side, agrees with SYM perturbation 
theory up to two loops. We proved the agreement for the general solution in the SU{2) sector. The two-loop 
agreement was also observed (for BMN states) at the level of 1/L corrections |SH1I^) which are described by 
quantum corrections in the sigma-model. In view of the discrepancies found in I42j for elliptic solutions at three 
loops, these results should be taken cautiously. Definitely more work is needed to understand if string solitons 
produce all-order non-perturbative results or the agreement is restricted to the first few orders of perturbation 
theory. We hope that the approach based on integrability is the right framework to address this question. 

It would be also interesting to generalize the approach developed in this paper to all scalar operators or 
even to the full set of single-trace operators in the SYM theory [211, for which the one-loop Bethe equations are 
known. It is natural to guess that the solution should be expressed in terms of differentials on more complicated, 
non-hyperelliptic Riemann surfaces. Again, we expect more complicated analytic structures to arise there, in 
agreement with the integrable structure of the full AdS^ x sigma models which were only briefly discussed 
in this paper. We hope to return to these problems elsewhere. 
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Note added: After we originally submitted this paper, we received |S2] where the transformation that maps 
Backlund charges in the sigma-model to those in the spin system was extended to two loops. We later received 
[5^ , where it was shown that the effective action for the Heisenberg chain can be mapped onto the string action 
to at least two-loop order. 
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A General finite-gap solutions of the sigma-models 

An alternative and more generic way to construct the finite-gap solutions for relativistic two-dimensional sigma- 
models, satisfying the Virasoro constraints, was proposed by Krichever in '84'. Here we review this approach 
and modify it in order to satisfy the Virasoro constraints (|4.3I) (with the non- vanishing r.h.s.). Moreover, some 
of the formulas of this construction appear to be useful for the examples we consider in sect. 15.21 and should be 
even more useful for the generalizations of our results beyond the SU (2) sector. 

The construction is based on an algebraic curve S with two fixed points P± endowed with a function 
X = x{P) which has simple poles at P = Qi, . . . , Qm where M is the complex "dimension" of sigma-model, i.e. 
the number of real co-ordinates in (|4.1(l is 2M. For the SU{2) case and sphere one should take Af — 2. 

These generic solutions are expressed in terms of the Riemann thcta-functions of F, which is the double-cover 
of S with the branch points P± , or in terms of the Prym theta- functions on S itself. The basic idea is that the 
sigma-model co-ordinates are determined by the values of the Baker- Akhiezer functions T{P,a±), satisfying 

{d+d-+u)T{P,a±)=Q, Per (A.l) 

where u is the Lagrange multiplyer of the sigma-model (|4.1|) , which ensures equality H4.2II . The solution for the 
sigma-model co-ordinates is given by 

Zi{(J±) = riT{qi,a±) Zi{a±) ^ riT{qi,a±) = riT{qi,a±), /=1,...,M (A. 2) 

where cr± = t ±a are relativistic co-ordinates on string world-sheet and Zi — X\ + 1X2, Zi = Xi — 1X2, etc are 
complexified target-space co-ordinates of the sigma-model. In l|A.2|l (q/, qi) are pairs of points on F conjugated 
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by the involution /, such that / : F i-^ F, = 1, and 

F/Z^E, f:P±^P±, 
f:{qi,qi)^iqi.qi). I^1,...,M (A.3) 
TTr-^s : (g/,?/) Q/, / = 1,...,M 

i.e. E is factor of F under the action of /, the points P_|_ and P_ are mapped by / to themselves, while the 
points qi and qj (whose projections on E are poles {Qi} of the function x on E are mapped by / to each other. 

The solution T(P, (t±) to (|A.1|I is the Baker- Akhiezer function on F with essential singularities at P± 
T(P, cr±) cx exp {n+{P)a+ + 9 {A{P) + V+a+ + U_ct_) - 

(A.4) 



p^Tp^ exp(fc±a±) {l + Y.ifia±)k^ 



where one can use the standard formulas (see, e.g. (2) ) with the inverse local co-ordinates k±{P±) — oo, the 
Abel map of F into its Jacobian A{P) — du), and U± — §-q dU,± are B-periods of the second-kind Abelian 
differentials on F 

dVL± ^ ±dfc± (1 + 0(fc^2)) , j dn± = (A.5) 

with the only second-order pole at P± respectively. 

The normalization factors in ljA.2|) are determined by 

r]^^^ , /=1,...,M A.6 

— x+ 

where dU, is the third-kind Abelian differential on E with the simple poles at P±, x± — x{P±) — resp^xdfl and 
normalizations (jA.Sp satisfy J^i rj = l due to ^ res (xdil) ~ 0. 

The Virasoro constraints (|4.3|l are ensured by 



(A.7) 



J2 \d±Zi\^ = ^rcsQ, ( d±Td±Tdn \ = resp^ ( d±Td±Td^ = 

J J \X— — Xj^ J \x+ — X- J 

f X — X± 2 \ 2 

= ztzresfe I =cxD fc+ ~; — I — K 

\x+ — X- K± J 

which is guaranteed after we require at P± 

xiP) ^r.^^ + (A.8) 

instead of used in "HT dx\p^ = which would be equivalent to k = in H4.3|l . For the sigma- model on S^, 
considered in sect. 14.21 we have M — 2 and the role of E is played by hyperelliptic curve H3.3|l . and the function 
X can be taken to be the hyperelliptic co-ordinate. 
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The double covers with only two branch points have relatively simple structure (see sect. 5 of for details). 
The period matrix of F is of the size 2g x 2g, where = g and gr = 2g, and has the form 

T f B + T B-T \ 
Tr = T:\ (A.9) 
2 \^ B-T B + T J ^ ' 

where T = is g x g period matrix of S and B is g x g Prym matrix of E. In the "non-relativistic" limit 
P± — > P, the cut on F between P+ and P_ shrinks to a single point, and F turns essentially into two copies of 
E. In this limit [HS| B = T + 0{a), where a is some parameter, measuring distance between P+ and P_, so the 
Prym theta- functions turn into the ordinary Riemann theta- functions on S. 

In the simplest case when E is a rational curve there are no theta functions and the Baker- Akhiezer functions 
are simply 

T = exp (fc_|-(7+ + fc_<T_) , T = exp (— fc+(7+ — fc_cr_) , (A. 10) 

where the local co-ordinates on sphere 

k+ = fc- - kVW, (A.ll) 



are expressed through a uniformizing function W, such that W{P+) — and VF(P_) = oo. Then the Abelian 
differential of the third kind is dSl — and generic form of the function x with two poles is 

^=E^^F^+^o (A.12) 

7=1,2 



(A.13) 



To get the desired form of expansion at zero and infinity ljA.8|) one should impose the constraint 

Wr \ ^ Wr I 

and with these data one gets from HA.6() 



1=1,2 ^ 



w+ - W- 



2 _ ^+~^ 



(A.14) 



r 



The function 1A.12II can be then written as 



If we parameterize sphere, as in (|5 



13 



y^^x +bx + c (A. 16) 



'Compare to sect. 15.21 we use here the coefficients b = B/A and c = C/A, (cf. with 15.411 ') and we recall that a - 
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the uniformizing function W can be chosen as 

W X — a-\-y — y+_z — 
Wo X + a + y — z — 

with yj^ — y{a), y_ — y{-~a), z± — y± ± a and a useful co-ordinate z is defined by 

z^ + bz + c 



2 

Z'' — C 



z ^ X + y, X = 



b + 2z' 



y = 



b + 2z 



The constant Wq is convenient to choose as 



and the poles in (|A.12|I . (|A.15|I to put into 



Wo' 



2 _ y+(6 + 2z_) 



y-{b + 2z+) 



+ 2z- 



(A.17) 



(A.18) 



(A.19) 



(A.20) 



corresponding to two infinities x — oo ljA.16p . 

The same function ljA.17|) with zero at a; = a and pole at a; = —a on the "lower" or "unphysical" sheet 
(corresponding in our notations to y = +\/ x'^ + fex + c) can be presented as 

W b + z+ + Z- z^ ~ z^ y + y^ 

or 



b + 2z- b + 2z- x + a 
Wo b + z++ z_ z+- z_ y + y+ 



W 

The resolvent (|5.41() can be presented as 

1 



b + 2z4 



b + 2z+ X — a 



(A.21) 
(A.22) 



G 



(6 + 2z+)(6 + 2z_) 



4(z+-z_) 
1/1 



y+y- 



4 \ X — a 



1 + ^] 



1 



y+ J x + a 



W 



y 



w 



i + — I I - 

y- 



_ iy+ - y-){b + z+ + z^) 
{y+ - y-)ib + z+ + Z-) 



(A.23) 



The normalization coefficients ()A.14|I are determined by (jA.Gp with x being "hyperelliptic" co-ordinate on (jA.16|) 

. y+-y- 1 , ^ , ^ T 

+ jj^ + 0{a), I = + 

(A.24) 



1 

2a 



dW 



2 4a 2 4Vc 

1 _ y+ - y- _ 1 b_ 

' 2 A^c 



4a 



where / = -!-,— and oo/ correspond to two points with a: = oo on ljA.16|l . From ljA.17(l one can easily find that 

+ (^{W ) ^ 2{z+ - Z-). — — — --^-f... 

^ ' ^ ^\l (b + 2z+)(b + 2z) kl 

^ , (A.25) 



a; -I- a 



x^a b + 2z+ Wo 
2y_{z+-z^) Wo 



+ 0(^-2) = _2{z+-z_) 



y+y- 



b + 2z_ W ' ' {b + 2z+){b+2z^)k'i 

i.e. the "hyperelliptic" co-ordinate x (the normalization here is inessential) indeed satisfies (|A.8p 
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B Elliptic integrals 



We use the following conventions for elliptic integrals: 



^(1 -a:2)(l -ra;2) ' 



(B.l) 



E{r) ^ / dx 



1 /I 2 

1 — rx^ 







Here we list some of their properties j8tij which are used in the main text. 

Under the modular transformation r — > 1/r elliptic integrals transform as follows 

K{l/r) = y^{K{r) -iK{l-r)), 
K{l-l/r) = ^K{l-r), 
E{r) + iE{l - r) - (1 ~ r)K(r) - irK(l - r) (b.2) 

E(i-i/.-) = Sl_l). 

The elliptic integrals of direct and dual moduli satisfy the Legendre relation: 

K{r)E{\ - r) - K{r)K{\ - r) + E{r)K{\ - r) = | . (B.3) 

C Reality condition for elliptic solutions 

The end-points of the cuts for the elliptic solution of Bethe equations must be complex conjugate of one another. 
This is not always the case for solutions of H5.77|l . H5.78I) . The reality condition thus puts a constraint on eligible 
m and n. We know that the solutions are real for m = or ji = (cases (i) and (ii) in sect. I5.4|l . since they 
reduce to the solutions found in 32 . Here we analyze solutions with arbitrary n and m and assume that to 7^ 0. 
We will analyze the reality condition in the limit of a 0. Then r is exponentially close to one and we can 
expand elliptic integrals as follows: 



Kir) = £ + 4(e - 1) e -2-^ + O (e^ e -^') , 

(C.l) 



K(\-r) = | + 2^e-2= + 0(e-4^) 



E{r) = 1 + 4(2£ - 1) e + O (e^ e -^^) 
E{\~r) = ^-27re-2^ + 0(e-*"), 

where r = l — IGe^^*^. Plugging these expansions into H5.77|l . we find 

mn 1 



^ 2to 2a 



l + 4e— +0 



(C.2) 
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and, from H5.78|l : 



Xl = 



X2 



la 
2m 
ia 
2m 



l-4e- 



1 + 4e- 



o(.-*) 



There are only two ways to satisfy the reahty condition. We can either put n 
which corresponds to case (i) in sect. 15.41 or we can take m — 2n, then xi ; 
(iii). 



(C.3) 



0, then Xl = — xi and X2 = — X2, 
— 5c2, which corresponds to case 
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